The N-chain Hubbard Model in Weak Coupling 
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We present a systematic weak-coupling renormalization group (RG) technique for studying a 
collection of N coupled one- dimensional interacting electron systems, focusing on the example of 
iV-leg Hubbard ladders. For N — 2,3, we recover previous results, and find that also more generally 
broad regions of the phase space of these models are unstable to pairing, usually with approximate d- 
wave symmetry. We show how these instabilities can be understood in terms of a fairly conventional 
"gap" function A at the discretized Fermi surface, and describe how this function is calculated. The 
dimensional crossovers as TV — > oo and as many such ladders are weakly coupled together are also 
discussed. 



I. INTRODUCTION 



Over the past few years, considerable interest has fo- 
cused on systems of coupled chain conductors. Early 
theoretical studies of Heisenberg ladders (appropriate for 
the strongly interacting, non-itinerant half-filled limit) 
revealed an interesting odd/even effect .EJ~El If the num- 
ber of legs of the ladder, N, is even, the system is ex- 
pected to be a spin-liquid with a singlet ground state 
and a (spin) gap to the lowest lying excitations carry- 
ing angular momentum. For odd N, the ground state 
has quasi-long-range antiferromagnetic order and a set of 
gapless spin-wave excitations, which puts it in the uni- 
versality class of the single spin-I/2 Heisenberg chain. 
Recent progress in the experimental preparation of rela- 
tively isolated spin ladders has begun to probe some of 
this rich physics and appears to have verified these ex- 
pectations for N = 2, 3.Erkl 

The behavior of doped ladders, i.e. those with itiner- 
ant charge carriers, is much richer. Particular theoretical 
attention hasj-been paid to the case N = 2, the two- 
leg ladder E'fJE3TL3 Early motivation stemmed from the 
possibility of realizing a cqnjCrete example of resonating 
valence bond (RVB) ideas .Oa According to this line of 
thought, since the two-leg Heisenberg ladder is a spin liq- 
uid, the doped carriers experience a short-range attrac- 
tive interaction, leading to pairing and the persistence of 
the spin gap. Such behavior is pdeed observed ir 
lations of two-chain HubbarcO^Eil and t-J models) 
which push the current computational limits of numerical 
methods working directly at zero temperature. Subse- 
quent work by numerous authors has since demonstrated 
the existence of such a spin gap phase for low dopings.by. 
controlled analytical methods in weak coupling .Ejo'Ej 
This weak-coupling approach has the additional advan- 
tage that it provides a full picture of the phase diagram, 



even away from half-filling. 

In this paper, we extend-thia analysis to more general 
TV-chain Hubbard models.LJCJ Such an extension is use- 
ful in two respects. First, it allows a determination of the 
phase diagram for any small value of N, thereby elucidat- 
ing the physics of doped spin liquids, the even/odd effect, 
and geometrical frustration. Furthermore, our equations 
allow a complete interpolation between one and two di- 
mensions (along a particular path in parameter space 
- see belpw). An understanding of such a dimensional 
crossovercil is a crucial first step in interpreting experi- 
ments in quasi-one-dimensional conductors. 

To determine the behavior in the weak interaction 
limit, we employ a generalization of the renormaliza- 
tion group (RG) developed in Ref. [f^ (the extension to 
the particular case N = 3 has already by studied by 
ArrigonScIl) . This provides a systematic basis for treat- 
ing the logarithmic divergences arising in a naive pertur- 
bative analysis. Coupled with the technique of bosoniza- 
tion, the primary output of the RG is a "gap" function 
A, describing pairing and the relative phase among the 
various spin and charge modes in the system. In the 
limit of large N, A becomes identical to the gap function 
defined in the conventional BCS theory of superconduc- 
tivity, and one may thereby connect our results directly 
with higher dimensional analogues. 

The RG also determines the zero temperature behav- 
ior as the chain length is taken to infinity. Because such 
a system is, for any finite N, still one-dimensional, it 
cannot sustain true off-diagonal long-range order, but 
is instead a generalized Luttinger liquid. The particu- 
lar Luttinger liquid phase, within a general classification 
scheme developed in Ref. [D] also follows from the gap 
function A. We will use this notation, in which a phase 
with m gapless charge and n gapless spin modes is de- 
noted CmSn, in what follows. 
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The results of our calculations for positive U Hubbard 
chains in the phase diagrams are summarized in Figs. [6- 
10]. We emphasize that the phase diagrams are valid for 
arbitrary filling n and transverse hopping t± except at 
some specific lines (see below). First note the prolifer- 
ation of phases as N increases from 2 to 4. We believe 
that this complexity persists even in the N — > oo limit 
(but see below) . The crossover to two dimensions is thus 
highly nontrivial. Secondly, despite the repulsive inter- 
actions, the majority of phases exhibit some degree of 
reduction of gapless spin modes, i.e. pairing. The sym- 
metry of the pair wavefunction is in most cases consistent 
with a d-wave form. Unlike the two-chain case, however, 
as N is increased, gapless spin modes exist due to the 
presence of nodes in the pair wavefunction. Under differ- 
ent circumstances, as can be seen from Figs. [6-10], both 
d x 2_ y 2 and d xy states appear. 

A number of special results are obtained for particu- 
lar small values of N. In the case N — 3, the differ- 
ence between open and periodic boundary conditions is 
quite pronounced, due to the strong role of frustration 
(and consequent absence of particle/hole symmetr jA, m 
the periodic-case. As found previously by Arrigonicj'tll 
and SchulzJ^j this Hubbard "prism" exhibits a spin gap 
at half filling, which persists over a range of both parti- 
cle and hole doping. An especially surprising effect oc- 
curs for N = 4 with periodic boundary conditions. In 
certain regions of the phase diagram, (singlet) Cooper 
pairs condense, not into the zero center-of-mass momen- 
tum state, but rather into the m = ±2 center-of-mass 
(quasi-) angular momentum states around the four-chain 
cylinder. We therefore call this a Cylindrically Extended 
(CEX) d-wave phase. Preliminary indications of the 
CEX d-wavfi-ehase have been found in recent numerical 
calculations SB 

Although detailed phase diagrams such as these have 
only been obtained for N — 2,3,4, our RG equations 
are valid for arbitrary N. They can be easily integrated 
numerically to any desired accuracy to determine most 
features of the weak-coupling phase diagram for any N. 
In the limit N -*-po, several connections can be made 
with other work.Ej This limit may be taken in several 
ways. For the simplest form of our flow equations to re- 
main valid, the interaction strength U must be scaled 
logarithmically to zero as N — > oo. Strictly speaking, 
then, these RG equations do not describe truly two- 
dimensional systems with finite, non-singular, interac- 
tions. The only logarithmic reduction of the domain of 
validity with increasing N suggests, however, that the 
two-dimensional limit may nevertheless be well approx- 
imated in this scheme. We present arguments that this 
is indeed the case. Firstly, in the large N limit, our RG 
equations reduce, up to an overall normalization of the 
interaction strength, to those of Shankar,Ej derived di- 
rectly in two dimensions. Second, an extended set of RG 
equations incorporating additional interactions, which we 
argue captures completely the two-dimensional limit for 
small non-zero U , can be shown to be equivalent to the 



previous ones for interactions with non-singular momen- 
tum dependence at the Fermi surface. 

Based on these analyses, we expect our RG equa- 
tions to contain a complete description of the dimen- 
sional crossover in the weak-interaction limit. In this 
limit, explicit analysis of the 1/iV correction terms show 
that all the pairing instabilities occur only at very low 
temperatures, T C (N) ~ Ajv ~ e^^. Feedback of the 
forward-scattering interactions into the Cooper channel, 
responsible for the pairing instabilities in the smaller lad- 
der systems at weak coupling, is therefore insufficient to 
produce superconductivity in the two-dimensional limit. 
We conclude that strong and/or nearly-nested interac- 
tions are necessary to explain superconductivity in 2d 
repulsive Hubbard systems. 

The remainder of the paper is organized as follows: 
In Sec. II, we introduce the N-chain Hubbard model, its 
weak-coupling limit, and a compact current-algebra no- 
tation for the allowed interactions. In Sec. Ill, the RG 
equations are derived using the operator product expan- 
sion to one loop order, and the numerical method used to 
study these equations is explained. In Sec. IV, we show 
how the results of these numerics can be understood using 
bosonization, introducing the gap function A in this con- 
text. These techniques are applied in Sec. V to determine 
the phase diagrams of 3- and 4-chain systems. Our analy- 
sis of the 2d limit is given in Sec. VI, and implications for 
numerics and experiments are discussed in Sec. VII. Four 
appendixes give further details of current algebra meth- 
ods, RG equations for umklapp couplings, initial values 
of coupling constants for the Hubbard models, and Klein 
factors needed for the bosonization calculations. 



II. N-CHAIN HUBBARD MODEL 



The N-chain Hubbard model is described by the Hamil- 
tonian H = H + Hjj , 

H *= E {-t[d\ a {x + \)d la (x)+h.c] 



-t 1 \d\ +la {x)d ia {x) + h.c. 
H U = Y J U: dl^d^dl^d^x) 



(2.1) 
(2.2) 



where d^dl) is a fermion annihilation (creation) opera- 
tor on chain i (i = 1 . . . N), and a =f, J. is a spin index. 
The parameters t and t± are hopping amplitudes along 
and between the chains, and U is an on-site Hubbard 
interaction. 

We begin by diagonalizing the quadratic part of the 
Hamiltonian, Ho, as appropriate in the weak-coupling 
limit, U <C t, tx- This is accomplished by transforming 
to new fermion fields tpi, where 
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E 



(2.3) 



The transformation matrix S depends upon the bound- 
ary conditions in the transverse (y) direction. For peri- 
odic boundary conditions (PBCs), the eigenfunctions are 
plane waves, and 



2wi 

S^^xl-eM^jm), (PBC) 



(2.4) 



while for open boundary conditions (OBCs), the trans- 
verse eigenfunctions are standing waves, 



5jm = VivTT sin( ivTT jm) ' (0BC) (2 ' 5) 

This brings the Hamiltonian into diagonal form in mo- 
mentum space: 



^° = E /^f^(P)4,0#ia(p), 



The single-particle energy is 



e i(p) — —2tcosp— 2t±cos(k y i). 



(2.6) 



(2.7) 



A difference in the spectra between OBCs and PBCs 
arises due to a difference in the set of allowed transverse 
momenta. These are 



kyl ~ N + 1 
2tt 



i, i = 1,2, ...,7V, for OBCs; 



(2.8) 



*v<=j7«> * = 0,±l,...,(±)[y] ) for PBCs, (2.9) 

where [a;] means the largest integer less than x. In the 
case of PBC, the momenta for k y i = ±ir are equivalent 
(i.e. differ by 2tt) for N ev en. For this reason, we have 
enclosed the final ± in Eq. 2.9 in parenthesis, which in- 
dicates that for N even, only one of these should be in- 
cluded for the proper counting of modes. 

Eq. 2.7 defines N bands, which, in weak-coupling, are 



filled up to the chemical potential (Fermi energy) /x. For 
those bands which are partially filled, this defines a set 
of Fermi points {kp i\ via 



(2-10) 



The chemical potential is fixed in terms of the physical 
density n (measured as a particle number per site) by the 
implicit condition 



(2.11) 



Z = Tre-> 3H = J [#][#]' 



(2.12) 



where (3 = {k^T) -1 is the inverse temperature. Unless 
explicitly stated otherwise, all calculations in this paper 
are performed at zero temperature, i.e. (3 = oo. The 
(imaginary time) action S is 



S = 



dr 



]ip ia (x)d T ^ ia (x) + H 



(2.13) 



and ip and ip are Grassman fields. 

In this formulation, it is straightforward to focus on 
the low-energy properties of the system. This is accom- 
plished by integrating out all Grassman variables cor- 
responding to fermionic operators creating excitations 
with substantial gaps. In particular, we first integrate 
out completely all ipia and ip ia corresponding to com- 
pletely filled or empty bands, for which all excitations 
are separated by a finite energy from the chemical po- 
tential. When interactions are included, this gives rise to 
modifications of the remaining effective action of 0(U 2 /t, 
U 2 /t±), negligible relative to the bare 0(U) couplings 
for U <C t. We denote the number of remaining par- 
tially filled bands by Nf < N. Secondly, we also inte- 
grate most of the longitudinal momentum modes in the 
partially filled bands, leaving only those in a width 2A 
around each Fermi point kpi (we will return to fix A 
later). That is, we integrate out ipi a {p) and 4>ia(p)i P ro ~ 
vided \p — kpi\ > A and \p + kpi\ > A. This again leads 
to a renormalization of the "bare" couplings in the re- 
maining effective action, this time with the additional 
logarithmic factor 



U 



1 + const, x — \n(k F /A) 



(2.14) 



This second step of integration is perturbatively con- 
trolled and makes negligible modification to the bare cou- 
plings, provided 



t 



ln(fc F /A)' 



(2.15) 



Assuming Eq. 2.15| is satisfied, the remaining fields 
have longitudinal momentum in only a narrow shell near 
the Fermi points. Within each shell, we can define chiral 
(right and left moving) fcrmions as 



for OBCs; (2.16) 



for PBCs. 



(2.17) 



We now turn to the treatment of interactions. It is use- 
ful to introduce a functional integral formulation. Corre- 
lation functions arc calculated as averages with respect to 
a "Boltzmann weight" e~ s , whose (Grassman) integral 
is the partition function (for example, see Ref. 



Here we have introduced the notation i = —i, which we 
will continue to use throughout the remainder of the pa- 
per. With this definition, ^Rii^pLi have opposite mo- 
menta in the PBC case (where transverse momentum is 
a good quantum number). The fields ipRi^Li may be 
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thought of as "slowly varying", due to their restricted 
range of momenta. 



For small A, the dispersion may be linearized within 
each momentum shell. The effective Hamiltonian is then 



^RiJ^Ria-AiJ^Lia > ( 2 - 18 ) 



where 



2t sin kpi. 



As it stands, the problem is formulated as an Nf- 
channel interacting Id Fermion system. It will sometimes 
be useful, however, to view the system instead as a finite- 
width strip of a two-dimensional Hubbard model. To 
translate between the two pictures, we recognize that in 
a hnite-size system, only a discrete set of transverse mo- 
menta k yi are allowed. One may think of these momenta 
as "cutting" the 2d Fermi surface, the intersections be- 
ing the Id Fermi points as shown in Figs. [1,2]. This gives 
an intuitive connection to more familiar two-dimensional 
physics, and also helps in identifying the possible four- 
fcrmion interactions. One caveat that should be kept 
in mind, however, is that for OBCs, the standing-wave 
transverse eigenfunctions are linear combinations of mo- 
menta ikyi, so that a single pair of Id Fermi points cor- 
responds in this case loosely to four points on the 2d 
Fermi surface. 
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FIG. 1: Band structure of the 4 chain Hubbard model 
with PBCs. Mapping onto the 2d BZ is shown at the 
left-hand side. The configuration of chains in real space 
is shown in the upper right. 
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FIG. 2: Band structure of 4 chain Hubbard model with 
OBCs. Since the actual transverse eigenstates with 
OBCs are standing wave superpositions containing both 
iky, we have indicated each Id Fermi point by two points 
(one closed and one open circle above) on the BZ. The 
spatial configuration of chains is shown in the upper 
right. 

We now try to write down all possible 4-point interac- 
tion terms allowed by symmetry. In addition to the U(l) 
x SU(2) symmetry corresponding to charge and spin con- 
servation, these terms must be preserved by charge conju- 
gation, time reversal, parity, and spatial translation op- 
erations. The most general particle-conserving 4-point 
vertex has the form 



Hi, 



A], 

Il^f E Sn(Q)V[{P a ,i a ,k a }} 



2tt ^ 

PaAa 



^ 1 u( fc l)V'k» 2 ( fc 2)V'P3 l 3(fc3)V'p 4U (fc4), (2.19) 

where Pi = ±1 «-> (R/L), and spin indices are left im- 
plicit. The fermion fields ipp a i a (k a ) are Fourier trans- 
forms of the slowly varying chiral fields defined in 
Eqs. 2.16-2.17. The strengths of the couplings are de- 
noted by T^[{P a , i a , k a }]. The total momentum transfer 
Q is 



Qx = -PlkFl - P2^F2 + P?.k F3 

—k\ - k 2 + k 3 + fc 4 , 



P A k 



2ir 

Qy = (-Pih 



Piii + Pzh 



(2.20) 

P 4 H) (PBC). (2.21) 



Note that Q y only appears for PBCs, since for OBCs, 
transverse momentum is not a good quantum number. 
Momentum conservation is implemented by the lattice 
delta function 



6u(Q) 



J2 nx 5{Q X - 2nn x ) (OBCs) 
En*,nJ {2 HQ-Mn x ,n y )) (PBCs) 



(2.22) 



where n x ,n y are integers. Vertices with non-zero n are 
called umklapp interactions. 
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The dependence of the vertex function on k a is an- 
alyzed via the Taylor expansion (we will assume it is 
differentiable) 



BV 

V[{P a , l a , k a }\ = V[{P a , i a , 0}] + V k a — + 0(k 2 ). 

^ Oka 



2A k y 



2A k y 



(2.23) 



We will see that, while the leading term is marginal in 
the RG sense, all the higher derivative corrections are 
in fact irrelevant, and can be neglected to the (leading) 
order of accuracy desired here. 



A. Interactions for OBCs 



Having obtained a general expression (Eq. 2.19) em- 



bodying the constraints on allowed vertices, we now turn 
to the classification of the solutions of these constraints 
in the particular cases of interest. We will do this first 
for the case of OBCs, proceeding in two steps. First, 
we locate the possible combinations of the band indices 
{Pa, ia}, and second, we determine possible combinations 
of spin indices, which are implicit in Eq. 2.19, by SU(2) 
symmetry. Constraints from other symmetries are also 
discussed. 

For OBCs, only momenta in the k x direction is con- 
served. We will assume A is sufficiently small that the 
internal momenta ki may be neglected in Eq. 2.2C. The 



condition for validity of this assumption will be derived 
at the end of the section. 

Since kpi < 7r, n x can take values, 0,±1,±2. For 
n x = ±2, all the fermi momenta kpi = tt. This means 
that these bands must be completely filled, and, follow- 
ing the reasoning described earlier, do not survive in the 
low-energy theory. For n x = ±1, momentum balance is 
possible in partially filled bands. However, at generic fill- 
ings, the fermi momenta are incommensurate (in units of 
2ir), and cannot be made to sum up to ±2tt. More care- 
ful consideration shows that such interactions only exist 
on specific umklapp lines in the (n, £j_/i)-plane. In this 
paper, we will restrict ourselves to the study of generic 
fillings, for which these umklapp interactions in the k x 
direction may be neglected. 




(a) 

FIG. 3: Examples of Cooper scattering c 
forward scattering /y (part (b)). 



(b) 

y (part (a)) and 



The last kind of vertices with n x ~ conserve x- 
momenta exactly. They may be found by plotting the 
interactions on the 2d Brillouin zone (BZ). These ver- 
tices satisfy 



Qx(n,tj_) = -P\k F i - P 2 k F 2 + PskF3 



PikpA = 0. 

(2.24) 



For a generic Fermi surface, two familiar classes of inter- 
actions are always allowed. The first comprises forward 
scattering interactions, which satisfy 



(P 1 ,P 2 )= (P 3 ,Pi) 
(ii>ia) = (*3j*4)i forward scattering, 



(2.25) 



The second set is the Cooper (or backscattering) channel, 
defined by 



Pi= P2,Ps = Pr, 
i\ = i%i is = ii, Cooper scattering. 



(2.26) 



In Eqs. 2.25-2.26, and in the remainder of the paper, 
P = —P and [xi, x 2 ) — {x3, X4) indicates pairwise equal- 
ity, i.e. either x\ — X3, x 2 = X4, or x\ = X4, x 2 = x%. 
The two possible solutions for forward scattering actually 
describe the same vertices, up to a sign from the Fermion 
ordering. Referring to the 2d BZ (see Fig. [3]), one sees 
that forward scattering interactions conserve the particle 
number separately in each band, i.e. one electron is anni- 
hilated and created in each of the bands i\ and i 2 . In the 
Cooper scattering channel, however, a pair of electrons is 
annihilated in band i\ and then scattered into band 

There are, however, other vertices at specific fillings. 
These vertic es correspond to special non-trivial solutions 
of Eq. [2.24 . Such solutions exist only on specific lines 
in the (n,t±/t) plane. Because these lines form a set 
of measure zero in the full phase space, the correspond- 
ing T^[{P a ,*a}] will be denoted "minor" vertices. Like 
the umklapp interactions, these minor vertices can be 
excluded at generic fillings. 

We have obtained the allowed vertices in momentum 
space. However, since the couplings are not momentum 
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dependent, they can equally well be written in terms of 
a local Hamiltonian density in coordinate space, i.e. 



B. Interactions for PBCs 



Hi 



F[{P a ,i a }]^ Pil ^ P ^ P ^ 



(2.27) 



where all the ip and ift are evaluated at the same space 
point. We have made the classification into forward and 
Cooper-scattering channels explicit by the change of no- 
tation V[{P a , i a , 0}] -> F[{P a ,ia}],C[{P a ,i a }], as appro- 
priate. In the case of forward scatte ring, as remarked 
earlier, the two solutions of Eq. 2.24 lead to a single F 
vertex. 

Each vertex obtained so far has several possible gen- 
eralizations once the spin indices are included. To make 
them explicit, we now introduce charge and spin cur- 
rents, in the scalar and vector representations of SU(2), 
respectively. These are 



When PBCs are imposed instead of OBCs, the sys- 
tem retains a finite set of discrete transverse translational 
symmetry operations. Correspondingly the transverse 
momentum k y (or more properly the exponential e lky ) 
is a good quantum number, and the allowed interactions 
are further constrained by the requirement Q y — 2-Kn y . 
Since kFi,k y i < it, n x ,n y can take the values, 0, ±1, 
±2. As explained in the last subsection, vertices with 
n x = ±2 can be ignored in the low-energy theory and 
those with n x = ±1 only live on specific umklapp lines 
and thus are not included. It follows that in the small 
U limit at generic fillings, it is sufficient to consider only 
vertices with n = (0, 0), (0, ±1), (0, ±2), which conserve 
the x-momenta exactly. 

The allowed vertices are found in two steps. First, we 
find all possible vertices which conserve x-momentum, 
and then we rule out some of them by conservation of 
y-momentum. To do the former, first note that under a 
(2.28) y-reflection k y — > — k y , which implies the Fermi momenta 



where <x denote Pauli matrices. These currents satisfy 
so-called Kac-Moody algebras, and this notation is there- 
fore often referred to as curr ent a lgebra. To regularize the 
composite operators in Eq. [2.28 , the currents are further 
defined to be normal ordered (although we do not indi- 
cate the normal ordering explicitly). Four-point vertices 
can be written down as products of two currents. Such 
Bilinear current-current interactions are SU(2) scalars (as 
appropriate for the Hamiltonian density) if and only if 
they are formed by couplin g two charge or two spin cur- 
rents. Each vertex in Eq. 2.27 therefore has two coun- 
terparts once spin is included. The subset of these which 
couple right and left movers is 



n 



(i) 



j 



-r p J R J L 4- r a T R 



rp jH TL , w JH . T 



(2.29) 



where and cy denote the forward and Cooper scat- 
tering amplitudes, respectively, between bands i and j. 
Summation on i,j is implied. 

Since fa , cu describe the same vertex, we choose the 
diagonal piece of the forward scattering amplitude to van- 
ish, i.e. fa — 0, to avoid double countings. Under charge 



conjugation Jy — * J 



ji, which implies 5y = 



Cji. Similarly, 

reflection symmetry (in x) implies — fji. While it 
is not obvious at this point, the choic e of signs for the 
scalar and vector couplings in Eq. 2.29 is such that they 
are all positive for repulsive on-site interactions. 

There are other interactions which are completely chi- 
ral, e.g. J R J R . As is well-known in conformal field the- 
ory, such purely chiral interactions do not renormalize 
or generate renormalization at leading order, and can be 
neglected in our weak-coupling analysis. Physically, they 
modify slightly the "velocities" of various charge and spin 
modes, which are already order one in the bare theory. 



satisfy kpt = k F j, where i = —i. This parity constraint, 
combined with conservation of a;-momentum alone allows 
vertices which satisfy 

(P 1 ,P 2 ) = {P 3 ,P i ), (ii,i a ) = (±i 3 ,±i4) (2.30) 



Pi=P 2) P 3 



P 4 



(2.31) 



Note that Eqs. ^30 
for OBCs (Eqs 
for PBCs. Physically. 



2.31 



differ from their counterparts 
2.25-2.26) by the extra choice of ± sign 
this arises because with PBCs, 
plane waves with momenta ±k y form two independent 
allowed transverse eigenfunctions, while only the single 
standing wave (superposition of the two) eigenfunctions 
satisfies OBCs. As before, additional vertices exist on 
special lines in the phase diagram, but will be ignored 
here. 

Since the ^-m omentum is also conserved, not all of the 
vertices in Eqs. 2.30[ -2.31 are allowed. Consider first the 
corresponding constraint for odd TV. We must evaluate 
Eq. ~ 



2.21, which can be rewritten as 



Piii - P 2 i 2 + Pais + PiU = n y N. 



(2.32) 



Since all the band indices satisfy |z| 
lutions with n y = ±2 do not exist. 



< (N — l)/2, so- 
Furthermore, after 
substitution of the partial solutions in Eqs. |2.30 



I2"32; i 



little algebra shows that the left-hand-side in Eq 
always even, so that no solutions exists for n y = ±1 ei- 
ther. Therefore, for the odd-chain systems, we need only 
consider the vertices which conserve momenta exactly, 
i.e. with n = (0,0). 



G 



Eq. |2.32| for n y = is satisfie d if and only if the 



sign is chosen in Eqs. pS(il -E> .3lL With this restriction, 
the allowed interactions are precisely the same forward 
and Cooper scattering ones that occur in the OB C ca se, 
and may therefore be described as before by Eq. 2.2£. 

The situation is more complicated for N even, because 
solutions of Eq. 2.32 exist with n y = ±1 and n y = ±2. 
The latter do not actually introduce any addition com- 
plications. This is because for n y — ±2, all the band 
indices must satisfy |«„| = N/2. The bands i a = ±N/2 
are, however, equivalent, since their y-momenta differ by 
2tt. Therefore, the band indices can instead be chosen 
equal, and then satisfy n y = 0. These n = (0, 2) ver- 
tices are thus included in the n = (0, 0) set which will be 
discussed later. 

The difficulty arises when n = (0, ±1), i.e. n y = ±1. 

- Pm - P 2 i 2 + P 3 i 3 + P 4 ii = ±N. (2.33) 



H 



(2) 



"''r l '' J jj 



-2a jR .T 



(2.36) 



iiji. Similarly, un- 
and R — > L, 



Impl icit in this notation is the constraint (see Eqs. 2.34- 
2.35 ) that the Uij are non-zero only for \i — j\ = -j- 

Since u^.iA describe the same interaction, we choose 

it 1 n 1 

m 2 - p , u 2 ? — to avoid double counting. Since under 
charge conjugation, Jy — > Jji, u, 
der a parity transformation, (i, j 
so Uij — ujj. 

Finally, of course, non-umklapp interactions also exist 
for even N. Just as for N odd, these n = (0,0) ver- 
tices are simply the forward and Cooper channel ones 
as for OBCs. Therefore, for even-chain systems with 
PBCs, the full set of allowed vertices comprises forward, 
Coope r, an d (Q y = ±27r) umklapp interactions, as given 
in Eq. O9] and Eq. |2.3€ 



In this case, choo sing t he — sign in Eqs. 2.30-2.31 leads 
to solutions of Eq. 2.32| . The first set of these is similar to 
the Cooper channel (see Fig. [4a]), with indices satisfying 



Pi — P2, P3 — Pi] 
h = 12,13 = U, 



k - h = ± ; 



N 



(2.34) 



The second set is similar to forward scattering (see 
Fig. [4b] ) , and has 



(Pi,P 2 ) = (P 3 ,P 4 ); 
(ii,i 2 ) = (h,h), H 

2A k y 



«2 = ± — ■ 



2A k y 



(2.35) 









JL I 


IjR 


kx 






, iR 








^ 




























y tl 


iR rlj 
































/ iR 

















kx 



(a) (b) 
FIG. 4: Examples of transverse umklapp scattering ujj 
(part (a)) and u?- (part (b)). As is clear from the figure, 
the momentum in the k y direction is not conserved. 

The two kinds of umklapp interactions in the k y direc- 
tion can also be written down as products of currents. 
Following the method we developed in last subsection, 
those y-umklapp interactions which couple the right and 
left movers can be described by, 



C. Constraints on momentum cutoff 

So far, we have determined all the interactions allowed 
by symmetry, assuming that the momentum cutoff A is 
"small enou gh" to neglect the — k\ — k 2 + k 3 + term 
in Eq. 2.2C. In this section, we make this requirement 
precise, and investigate how this begins to break down 
for larger A. 

In fact, the results in the last two subsections are 
strictly correct only for A = 0. With a finite cutoff, 
the picture is modified in two ways. Firstly, the specific 
lines on which the "minor" couplings exist are widened 
and occupy a finite area in the (n, t±/t) plane. Secondly, 
for sufficiently large A, new vertices (not included in the 
forward and Cooper scattering channels) can arise for 
generic fillings (i.e. throughout the (n,t±/t) plane). 

Consider first the minor vertices, which exist in the 
region of the phase diagram defined by 



f(n, tj_;ijkl)= -Pikpi - P 2 k Fj 



f(n,tj_;ijkl) 



< 4A. 



P3k F k + Pik F i, 

(2.37) 



For A = 0, as noted previously, solutions of this equation 
other than the Cooper and forward scattering ones exist 
only on isolated lines. Since f(n,t±) is a smooth func- 
tion of its arguments, there will be fi nite n eighborhoods 
around these lines which satisfy Eq. |2 . 37 . The widths 
Sn of these neighborhoods can be estimated by Taylor 
expanding / around the / = lines, i.e. 



dn 



A. 



(2.38) 



Since the derivatives of the fermi momenta, dkp/dn, are 
order one, the width of the line is approximately as large 
as the cutoff, 



Sn ~ - — . 

k F 



(2.39) 
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To know the fraction of the phase diagram influenced 
by these minor vertices, it is also necessary to determine 
the number (or equivalently, the density) of these regions. 
Consider first OBCs. Since each region grows adiabati- 
cally from a A = line, we can simply count the num- 
ber of solutions to f(n,t±;ijkl) — (at, say, fixed 
Roughly, the number of solutions may be estimated as 
follows. Picking a fixed n and t±, we choose three of the 
band indices, say i,j and k. Then / = determines a 
Fermi momentum kpi for the last band. Generally, how- 
ever, this momentum will not be one of the discrete set of 
Fermi momenta for this n and t± . Now begin varying say 
t±, keeping i, j, k and n fixed. As t± varies, so do kFi,kFj 
and fcpfe, and hence the required /cfz- As this happens, 
very soon kFi will pass through an allowed value, and we 
have found a solution. Given that, one may then vary 
both t± and n in order to keep / = for this particular 
ijkl, defining a curve in the (n,t±/t) plane. Since this 
can be repeated for each set of i,j,k, the total number 
of such curves is N minol ~ N 3 . For PBCs, roughly the 
same argument holds, except that k and I are related by 
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(a) (b) 
FIG. 5: New vertices for finite cutoff A. It is possible to 
shift the Cooper vertices (shown by dashed lines) for a 
finite cutoff and still maintain momentum conservation. 



The maximum shift S is estimated in Eq. 2.44 



Let us consider in detail the conditions under which Tt s 
conserves momentum in the cut-off theory (see Fi g. [5]) . 
Because both i and j have been shifted by 6 in Eq. |2.41 , 
conservation of Q y . The number of minor vertices for Q y is already conserved by design. Taking into ac- 

A^ 2 . cou nt Q x conservation by Taylor expanding (for 8 <C N) 

Eq. 2.37 gives the requirement 



PBCs thus reduces to Af„ 



For large N, the fraction of the phase diagram in which 
minor vertices contribute can be estimated simply by 
summing the widths of these lines. The resulting frac- 



tion /„ 



N„ 



T 5n is negligible (/ <C 1) provided 



A 1 

k^^N 3 
A 1 

kF~ < iV 2 " 



(OBC), 
(PBC). 



(2.40) 



dk-^-dk- y ^ 



— 5 < Aa. 

N 



(2.42) 



To obtain an order of magnitude estimate, we then ap- 
proximate the mean curvature between bands i,j by its 
typical value tj_/t, 



dh, 

d,k 



J\ 



d 2 k x 2tt 
tj_ 2n 

— — * — 7 • 

t N 



(2.43) 



The maximum allowed shift ^ max in band indices is there- 
fore 



For sufficiently large A, the one-dimensional bands as- 
sociated with adjacent points on the Fermi surface begin 
to overlap, and it becomes possible to form new inter- 
actions by substituting one for the other in the original 
forward and Cooper scattering channels. If A is large 
enough to allow this, such interactions exist generically, 
i.e. throughout the (n,t±/t) plane. As an example, con- 
sider the shifted terms (for PBCs) shown in Fig. [5]: 



i+SJ+S'Jij + fiji^Ji 1 



J L 

i,i+5°j,j+8- 



(2.41) 



t ^ N 2 A 
tx \i-j\k F 



(2.44) 



For (5 max < 1, only the unshiftcd Cooper channel interac- 
tion is allowed. Demanding this leads to the constraint 



A ,*]_., 1 , 



(2.45) 



Combining the constraint on the initial coupling 
strengths, Eq. 2.15, and that on the momentum cutoff, 



Eqs. ^40|, ^45T we find the reduced set of interactions 
in Eq. 2.29| and 2.36 is sufficient provided the initial cou- 
pling strengths satisfy 



Here S parameterizes the transverse momentum shift; for 
5 = 0, the vertices reduce to the familiar Cooper and 
forward scattering types. 



U 1 

J <C h777' 



(2.46) 



Here we have dropped the order one factors in front of the 
logarithm which are different for OBCs and PBCs. Since 
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the condition on U is only logarithmic in N, it is not a 
severe constraint on the initial values of the couplings for 
finite chains. Nevertheless, it is clear that the true two- 
dimensional limit is rather subtle. Indeed, Eq. 2.46 indi- 
cates a possible non-commutativity of the order of limits 
(at zero temperature) in taking U — ► and N — ► oo. 
We will return to the interesting and important issues 
involved in the 2d limit in Sec. VI. For the moment, we 
will restrict ourselves to finite N and discuss the corre- 
sponding weak-coupling behavior of such systems, under 
the conditions of Eq. 2.46. 



III. 



RENORMALIZATION GROUP FLOW 
EQUATIONS 



Ia 



Each 



using the compact notation _ ja < \ Xi _ X2 \ < b- 
integral is now separated into two parts: long-wavelength 
modes, \x\ — x%\ > ba, and short-wavelength modes, 
a < \xi — X2\ < ba, where b > 1 is the rescaling pa- 
rameter. This separation is convenient because in the 
latter integral, all fields are at nearby space points. This 
allows the use of the operator product expansion (OPE) 
to replace the products of such nearby operators by a 
series of local operators (for review, see Ref. [}?],[38]), i.e. 



Y\_ dr i dx l (Hint (ti,Xi )H in t (t 2 ,x 2 ))q 



i=l,2 



dxdT(8TLi 



(3.4) 



To analyze the behavior of the weakly interacting sys- 
tem, we employ the RG approach. In this section, we de- 
scribe the scheme used and present the resulting differen- 
tial RG flow equations (analogous flow equations, for par- 



ticular restricted cases were obtained in Refs.E3 andE3). 
Further details of the calculations can be found in ap- 
pendix A. The general approach of the RG is to progres- 
sively eliminate short- wavelength, high-energy degrees of 
freedom. To formulate this mode elimination, we first 
rewrite the partition function in terms of an average, 



Z= / [#][<#] 



-So-Si. 



= Z (e' 



(3.1) 



where Zq is the partition function without interactions, 
and angular brackets with the subscript denotes an av- 
erage with respect to the quadratic action So only. This 
form may be re-exponentiated using the cumulant expan- 
sion, 



1, 



) = CX P <{ (-Sint)o + ^[(Sf„ t >0 - (Sint)o] 



+o(s? nt : 



(3.2) 



Up to this point, we have systematically derived the 
low-energy Fermion model with a "momentum-shell" cut- 
off A. While the RG may be implemented directly with 
this model, it happens that the one-loop RG equations 
needed here are in fact independent of the cut-off scheme 
used. This independence arises from the dominance of 
logarithmically divergent terms at one-loop level, whose 
coefficients are insensitive to the particular form of cut-off 
used. We take advantage of this property here by adopt- 
ing instead a real-space cut-off a w 1/A. This distance 
then appears as an explicit cut-off in all x-integrals, e.g. 



poo 

>0= J\ dT i dx i (Hint{Tl,X 1 )Hint{T2,X 2 ))o, 

Ja i=l,2 



(3.3) 



The method to compute the effective interaction, d7ii n t, 
can be found in Appendix A. 

As shown in appendix A, the effective interaction 8TLi n t 
has the same form as the original Hamiltonian, and thus 
has the effect of renormalizing the bare couplings. The 
RG finishing with a rescaling step, which attempts to 
bring the theory as much as possible back to its origi- 
nal form. To restore the original value of the cut-off and 
maintain the original set of Fermi velocities requires the 
change of variables 



(3.5) 



, X , T 

_ x = V T = V 

^(x, t) = b 1/2 ip(x', t'), ip'(x, t) = b 1/2 i;(x', t'). (3.6) 



While this indeed preserves (at one-loop order) the form 
of 5*o, the interactions are of course changed. The sim- 
plest way to keep track of these changes is to perform the 
RG inEnitesimally, with the rescaling factor b = e dl . It- 
erating both steps of the RG then leads to differential RG 
flow equations for the coupling constants, as a function 
of length scale L(l) — e l . 

For OBCs, the allowed couplings are forward and 
Cooper scattering. The RG equations governing them 
are 



ff.-. 



i4 



3 

16 



(4) 2 , 



-(/, 



cr\2 



(3.7) 
(3.8) 



+(c p h p + —c? h? ) 



(3.9) 



(3.10) 







where /y = fij/n(vi + Vj) and the same for c y -. Also, 
we define /ly = 2/y + JyCji for convenience. The weight 



factor in the summation a,- • = ( t, '+ 1,fc )( t, 3+ t ' fc ) 



2v k (vi+Vj) 



is sym- 



metrical in i, j. The dots indicate logarithmic derivatives 
with respect to the length scale, i.e. / = df/dl. 

For PBCs, if the number of chains is od d, the all owed 
vertices are the same as in OBC, and Eqs. |3~7- 3.1C hold 
without modification. However, if the number of chains is 
even, the additional umklapp interactions in the k y direc- 
tion give rise to additional flow equations and extra terms 
in the forward and Cooper scattering equations above. 
Because these augmented RG equations are quite com- 
plicated and do not provide any obvious insights upon 
inspecti on, we h ave relegated them to appendix B. 

Eqs. [3~7| - |3.1C| must be supplemented by initial con- 
ditions to completely specify the problem. While it of 
course straightforward to obtain these initial values from 
the bare couplings, it does require a certain amount of al- 
gebra t o work o ut the effect of the unitary transformation 
in Eqs. 2.4-2.5 (see appendix C). Initially, then, all the 



couplings fij,Cij are 0(U/t) (but with specific ratios) 



IV. STRONG-COUPLING ANALYSIS 

In this section, we describe the application of the RG 
equations to the JV-chain Hubbard models with both 
OBCs and PBCs for small N. We will see that the RG 
flows almost always diverge, and will discuss the interpre- 
tation of such divergences, using as input certain general 
features of the numerical integrations. The instabilities 
encountered generally correspond to some degree of pair- 
ing, the notion of which we make precise through the 
study of various pair fields. 



A. Classification of couplings 

1. General scheme 

With the initial values in hand, the RG flow equa- 
tions can be integrated to investigate the physics of the 
weak-coupling limit. We have done this using Mathe- 
matica on a Sun Sparc-4 workstation. While the specific 
solution found depends upon the details of the model pa- 
rameters (e.g. N, t±/t, n), certain gross features of the 
behavior are generic. In particular, for almost all sets 
of initial conditions, the solutions of the RG equations 
are singular, and certain linear combinations of coupling 
constants diverge at some finite Id- Since the RG equa- 
tions were obtained perturbatively, they are not valid 
arbitrarily close to such an apparent divergence. To 
obtain sensible results, we instead cut off the RG flow 
at some specific length scale I* < Id, chosen so that 
U/t -C ma,x{ fij,Cij,u.ij} <C 1. At this cut-off length 
scale, the couplings may be classified into two groups. 
The first set includes those couplings which have become 



"large" but still weak, i.e. U/t <C gi <C 1, which we 
call marginally relevant. The remaining couplings do not 
grow under the RG, but remain 0(U/t) or smaller, and 
will be called marginal or marginally irrelevant respec- 
tively. At the length scale I*, the system thus exhibits a 
separation of energy scales, with the marginally relevant 
interactions much larger than the marginal or marginally 
irrelevant ones. The phase diagram of the system may 
then be determined simply by neglecting the latter in- 
teractions and studying the states determined by the 
marginally relevant couplings alone. 



2. Strict U -> 0+ limit 

In the truly asymptotic limit U ^ + , much of the 
classification of couplings can be accomplished analyti- 
cally. To do so, consider the formal solutions of the RG 
flow equations as functions of I and the Hubbard interac- 
tion U, which we will denote gi(l; U), where i is a compos- 
ite index labelling all the interactions. The perturbative 
RG is valid provided all the renormalized interactions 
are small, i.e. \gi(l; U)\ <C 1, which is certainly true ini- 
tially. To proceed, we need to assume something about 
the form of gi(l; U) near the divergent scale Id- Assuming 
a power-law first suggested in Ref. 30 (which can be ver- 
ified analytically for simple cases and numerically quite 
generally), dimensional analysis essentially requires 



9i(i;U) 



(la - Ulp* 



(4.1) 



where Iq = Uld and we have set t = 1 here and in the re- 
mainder of this section. The order one coefficients Gi and 
exp onen ts 7.; remain to be determined. Given the form of 
Eq. 4.1, it is clear that the formally divergent couplings 



only become much greater than U when I is very close to 
Id- We are thus actually interested in the asymptotic be- 
havior of Eqs. 3.7-3.10. Unlike the full integration of the 



RG flows, the less ambitious task of determining these 
asymptotics can in fact be accomplished analytically, as 
we now demonstrate. 

In order to fix the exponents 7j, it is necessary to use 
some input from numerics. In particular, in every case we 
have examined, the vector part of the forward scattering 
interactions, is always marginally irrelevant. From 



Eq. 3.8, this implies that the combination 2c; 



1 \ y / 

is small, i.e. < 0(U 2 ). Since we are keeping only those 
interactions which scale to values of order one, to this 
order of accuracy the renormalized couplings satisfy 



mi 



0, 



(Mi) 



(4.2) 
(4.3) 



Eqs. 4. 2 -4. 3 suffer corrections of 0(U), but may be treated 
as equalities in the following leading-order analysis. 

Eq.|4~2}|0 can be understood in a simple physi cal way. 
Some simple algebra demonstrates that if Eqs. 4.2-4.3 
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were replaced by exact equalities, these conditions would 
be preserved by the RG flow. Thus it is natural to suspect 
that under these conditions, the system has acquired an 
additional symmetry. Indeed, upon closer examination 
one finds that the constraint implies independent conser- 
vation of spin within each band. Although this is not an 
exact property of the Hubbard model, it is approximately 
satisfied due to the on-site nature of the interactions. For 
on-site interactions, Fermi statistics allow only a cou pling 
of oppositely oriented spins, which implies Eq. 4.3. Ap- 



parently the deviation from this symmetry caused by the 
nonzero initial values of is sufficiently small to allow 
the symmetry to be asymptotically restored at long dis- 
tances. 

Based on this observation, we will calculate the ex- 
ponents 7i in the U — > + limit, taking as an example 
the Cooper scattering vertices. Other exponents can be 
obtained by similar means. Th e RG equations for the 
Cooper couplin gs c f t in Eq. 3J3 can be re-written with 
the help of Eq. EO as 



The relative strengths of various couplings in the asymp- 
totic regime can then be determined (relatively) easily 
using these algebraic relations. Usually there is more 
than one solution for Gi. The identification of specific so- 
lutions with specific initial conditions which depend on, 
e.g. filling factor n, can only be found using numerical in- 
tegration of the full RG equations. We emphasize as well 
that these results rely upon on the strict U — > + limit. 
A priori, given the lack of other energy scales in the Hub- 
bard model, we would nevertheless expect such results to 
hold qualitatively provided U < 1 (for fixed small finite 
N). In fact, given the large number of coupling constants 
involved, apparently order one factors can conspire to 
render the limits of validity of this strict weak-coupling 
limit considerably smaller (e.g. U/t < 10 -5 in some re- 
gions of the phase diagram even for N = 2 - see section 
V!). For U/t < 1 but outside the strict U -s- 0+ limit, 
the algebraic relations do not hold, but we nevertheless 
expect quantitatively correct results from the numerical 
integration of the full RG flows. 



(4.4) 



Since all terms on the right-hand side are of the same sign 
(negative), they cannot cancel one another, and balanc- 
ing the two sides of the equation then gives the constraint 



where yff = j(c[ 



l = max{2 7 -,2 7 ^} 



Eq. 4.5 has two solutions. The first is 



ctr c 

7* < 7« 



1. 



(4.5) 



(4.6) 



B. Bosonization strategy 

Our strategy will be to bosonize relevant couplings in 
the renormalized Hamiltonian at the. scale I*. This is 
done using the Bosonization formulao'Lj 



(4.9) 



Here the chiral Boson fields 4>R/Lia obey the commuta- 
tion relations 



The second possibility is, 



(1+7D/2<1- 



(4.7) 



We thus conclude that all the exponents associated with 
Cooper couplings are bounded above by one. Similar 
considerations applied to the other RG equations imply 
that all the exponents are less than or equal to one. If 
one can probe arbitrarily near the divergent point Id (i.e. 
if U is put arbitrarily close to + ), any couplings with ex- 
ponents 7i = 1 will eventually outstrip any others with 
smaller exponents (even those with larger prefactors). In 
the strict t/ — > + limit, therefore, the relevant couplings 
are those with exponents 7j = 1. 

From Eq. 4.1, the relevant couplings with initial inter- 



action U thus satisfy 

9i(i;U) 



UGj 

ln — UV 



(4.8) 



for I near the cuto ff length scale 
into Eqs 



3.7-3.10 



. Substituting Eq. \L8 
the parameters U and Id cancel out, 
leaving a set of algebraic relations between the constants 
Gi. These alg ebraic equations are formally obtained from 
Eqs. 3.7-3.1C simply by replacing gi — > Gi and gi — ► Gi. 



[<t>R%a{x),<f> Rj p{y)\ = -[<j>Li a {x),<j)L fi{y)] 

1 

= -sgn(x - y)5ij5 a p, (4-10) 

1 

[4>Ricx{x),<t>Ljp{y)]= 4<Ma/3 (4.H) 

The rji a are Majorana (real) Fermions, known as Klein 
factors, introduced to preserve the proper anticommuta- 
tion relations between Fermion fields with differing band 
and spin indices. They obey 



(4.12) 



It is usually more convenient to trade the chiral Boson 
fields pairwise for a conventional Bosonic phase field <$> 
and its dual (displacement-like) field 9, defined by 



VRia + 4>Lia 
^Ria — <t>Lia- 



(4.13) 



They satisfy [(f>(x),8(y)] — —i®(y— x). Physically, the 
8(x) field describes the displacement of the electrons, 
while the dual field 4>{x) represents their phase. We can 
make a further canonical transformation to 
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(0,0)<„ = [(0,0)iT-(0,0k)]/V2. 



(4.16) 



(4.14) 



The p fields then describe charged singlet excitations, 
while the a fields describe neutral excitations carrying 
spin. Carrying through the change of variables in Eq. 4.9 
carefully, one finds that the non-interacting Hamiltonian 
Eq. 2.18 is equivalent to the Bosonic Euclidean action 



where 



X.T * 



[(d x cj> iv f + (d x 6 iu ) 2 ] + id x 6 iv d T ^ v , 



(4.15) 



p,a. 



Upon bosonizing with Eq. 4.9, the four-Fermion inter- 
actions are converted to linear combinations of gradient 
couplings and sinusoidal functions of the phases. The for- 
mer give rise to continuous shifts of the parameters of the 
low-energy description: modifications of mode velocities, 
charge stiffnesses, etc. We say that these shifts leave the 
system in the same phase. The sinusoidal interactions, 
by contrast, can potentially cause more drastic changes 
in the low-energy theory. They tend to "pin" their argu- 
ments (linear combinations of the Bosonic phases) to par- 
ticular values, modulo 2ir. We will treat the marginally 
relevant couplings in this way, by expanding the corre- 
sponding harmonic functions around their minima, re- 
garding the fluctuations around them as massive. This 
will effectively give gaps to some of the modes of the orig- 
inal non-interacting Fermion system, resulting in distinct 
new phases. 

To perform this procedure and determine the nature 
of the gaps arising from the interactions, then, it is suf- 
ficient to keep only those marginally relevant couplings 
which become sinusoids upon Bosonization. There are 
still quite a number of potential such terms, so it now 
pays to use some additional input from the numerics. 



C. Generic Instability 

For simplicity, we will first focus on situations in 
which the relevant couplings include only forward (/) 
and Cooper (c) scattering vertices. This is completely 
general for OBCs, but excludes certain regions of phase 
space in the case of PBCs with even N, for which the 
transverse Umklapp (it) couplings may become relevant. 
This alternate channel of instability will be returned to 
later in this section. 



1. Semiclassical Hamiltonian and analysis 

As a preliminary step in the analysis, we first rewrite 
the interactions in terms of the underlying Fermion fields. 
Using the SU(2) identity, 



the scalar and vector parts from Eq. 2.29 become 



H 



(i) 



E 



fij ^Ria ^Ria^LjP^LjP 



H4j + \ c 1j)^R la ^R 3a ^Li^Ljl3 



(4.17) 



Here we have dropped the irrelevant forward scattering 

j ij 

Following the strategy above, we must now Bosonize 
the system. In doing so, the interactions in Eq. 4.17 may 



be divided into two sets. The first consists of fr. and c^, 
which only contribute gradient terms after bosonization. 
To determine the phase of the system, therefore, we need 
include only the second set, which contains the Cooper 
couplings c£, cf , and cf.-. Keeping only these terms, 



and imposing the constraint in Eq. 4.3, the interaction 
Hamiltonian becomes 

^■int = 2 C ^ i ^ Ria ^ Lia ' l ^ Ria ^' Lia E 

^Ria^Lia^Rja^Lja + ^Ria^Lia^Lja^Rj^ > ( 4 - 18 ) 

where a = —a. The inter actio ns with a = (3 in the third 



and fourth terms in Eq. [4 . 1 7| cancel each other, giving 
rise to the aforementioned independent conservation of 
spin in channels i and j. This ordering of fermion fields 
is particularly convenient for Bosonization. 

As shown in detail in Appendix D, for this set of in- 
teractions the Klein factors can be repre sente d by the 
identity r/i = 1. Inserting Eq. 4.E in Eq. 4. IS , we then 
obtain 



(i) 



E 



c*cos(V8^M + E 4c ° 



U 



xcos(V4tt^ )cos(V27r0 i(T )cos(V27r6' ; ,- (7 ), (4.19) 

where <%f = ± ^)/V2. As pointed out by Schulzjl 
this explicit form after bosonization depends on the rep- 
resentation chosen for the Klein factors. In other words, 
we still have some "gauge" freedom left to shift the 
bosonic fields. However, the physical correlation func- 
tions, which include these Klein factors, are independent 
of the specific gauge choice. 



We next locate the minima of Eq. |4.19| . These can a 
priori be nontrivial, but turn out to be very simple in 
practice. In fact, we find that in all cases, the numeri- 
cally determined values of t he co efficients and are 



such that each term in Eq. [4. 19| can be minimized sep- 
arately. Most often this occurs because all the relevant 
couplings occur in channels connecting only two specific 
bands. We will focus on this special case now, in order 
to present a more detailed but containable exposition. 
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To proceed, let us denote the indices of the two 
strongly interacting bands by a and b. The numerical in- 
tegration demonstrates that, although they are initially 
positive (repulsive interactions), the diagonal Cooper- 
scattering spin vertices are driven negative under the RG. 
That is, 



Ci c bb < 0- 

We also obtain the sign 



c ab — ^ c ab ^ U ' 



(4.20) 



(4.21) 



Given these signs, a global minimum of Eq. 4.19 is 



2tt(0, 
2^ 
4W 



ba) — 
~P\ - 



LIT, 

77171", 

(I + 777 - 



2?7 + 1)7T, 



(4.22) 
(4.23) 
(4.24) 



where I, m, n are integers. Since all solutions give the 
same results of correlation functions, we will pick the 
I = 777 = 77 = solution for convenience. Fluctuations 
around this semiclassical solution are massive, as can be 
seen by the change of variables 



lao = 


Qacr) + S9 arT , 


(4.25) 


%* = 


8ba) + S8bcr, 


(4.26) 


P ab = 


€b) + s€b- 


(4.27) 



Expanding to quadratic order gives 

+\{mlf{SB Ur f + l(m p J 2 (Sr-) 2 , (4.28) 
up to a constant. The masses for spin and charge modes 



1/2 



m° b =2V2^(\4 b \+c: b ) 1/2 



C = V^« 6 ) 1/2 . 



(4.29) 
(4.30) 
(4.31) 



Comparison with the corresponding quad ratic frequency 
terms in the non- interacting action (Eq. 4.15) indicates 
that Eq. 4.28 describes three gapful modes, 6 aa , Q b(J , and 
. The first two terms in Eq. 4.28| suppress fluctuations 
in the spin-densities of the two channels, and correspond 
to spin gaps with magnitudes proportional to m" a and 
m b . The third term "locks" together the relative phase 
of the charge modes in the two channels. The remaining 
linear combination <^J~ and its conjugate 9 p b are not af- 
fected by Hj^ t , and continue to describe a gapless total 
charge mode. 



2. Pair fields 

The existence of a spin gap naturally suggests pairing 
of electrons with oppositely oriented spins. To investi- 
gate this notion further, it is instructive to consider a 
pairing operator 

i>(X,x)=iJ>Pi a (X + x/2)iPp, jl 3(X-x/2), (4.32) 

which annihilates two electrons with specified band in- 
dices and spin at particular positions. For compactness, 
we will omit explicit labeling of the pairing operator if? 
unless necessary to avoid ambiguity. In a true super- 
conductor, such pair fields condense, so that ^ 0. 
The Mcrmin- Wagner theorem prohibits such a continu- 
ous symmetry breaking in 1 + 1 dimensions, but \t can 
have power-law correlations (quasi-long-range order). Of 
physical interest are correlation functions of two or more 
of the various if? operators at different well-separated 
points (in X). For instance, 



C AB (X,x) = (& A (X,x)i? B (0,x)), 



(4.33) 



where the composite index A = {Paiac^a, P'aJaPa), and 
B is defined in the same way. Following the RG strat- 
egy, correlation functions such as Cab(%) are evaluated 
in several steps. First, we employ the perturbative RG, 
integrating out Fermion modes until the scale I*. Ignor- 
ing perturbative corrections from the mode integration, 
each pairing operator then picks up just the rescaling 
factor 



1 



*(X/A£,x/A£ ; r), 



(4.34) 



-v* 



where we have defined the coherence length £ = A~ 
At this point, the relevant couplings have become of or- 
der one, and may be safely bosonized. Carrying this out 
gives 



^(X,x) ~ — r)i a r) jf 3 exp 



iV^(<t>p ia [(X + x/2)/A£\ 
+<t>p, jft [{X-x/2)/Mi]) 



(4.35) 



The next step is to insert this in the desired correla- 
tion function, 



e.g. 



Ec 



4.33 , and integrate out the mas- 
sive modes using Eq. 4.28. Since the masses in Eq. 4.28| 
are order one, the correlation functions of the gapped 
phase variables will decay exponentially over (rescaled) 
distances of order one. This allows us to make the re- 
quirement of "well-separated" points more precise. Since 
the internal coordinates of the fields have themselves 
been rescaled by £, it follows that for \X\ 3> £, the mas- 
sive phase variables in two pair fields separated by a dis- 
tance | JSC | are exponentially decorrelated (provided the 
"internal" coordinates satisfy |x| <C \X\), and may be 
integrated out independently for each *f? . For instance, 
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& A (X,x)* B {0,x f )) k{& a (X,x)) (* b (0,x') 

\ I g \ I g \ I g 

(4.36) 

for \X\ ^> £, where the subscript g indicates an average 
over the gapped phase fields. 

It is thus sufficient to study the partially averaged pair- 
ing operators ($(X, x)) g . We will do this carefully for the 
case of the two bands a and b. The average is carried out 
with respect to the action S = So + S% , where 

So = [ ?[(W 2 + [dJsf] + id x 9 s d T <p s , (4.37) 
where s — aa, ba, abp—, as obtained from Eq. 4.15| , and 



Si = 



(4.38) 



From Eqs. 4.28 and 4.37, the small deviations 86 aa , 86b a , 
and 8<j) ab ~ are decoupled (from each other, but not from 
their conjugate fields) at the quadratic level, and may 
therefore be averaged independently. 

In carrying out this average, it is important to note the 
appearance of the conjugate fields 8cj) aa , 8<ftb<r, and 86 p ab . 
By the uncertainty principle, since [4>,6] = 0(1), these 
variables are wildly fluctuating. This implies that any 
complex exponential containing one of these fields will 
average to zero, unless it appears in the form of a "neu- 
tral" difference at nearby points, in which its fluctuation 
mean value is automatically subtracted. If such a sub- 
traction does occur, the average will decay exponentially 
with the separation of the subtracted fields. 



In fact, only four pairing operators satisfy this strong neutrality constraint, and are therefore non- vanishing. These 
are 



(* d+ (X,x)) g = (AZ)-^ndrKX + x/2)/AZWLdi[{X-x/2)/AS]) g> 
(4> d -(X,x)) g = {AZ)- l (1> m [(X + x/2)/A£]1,Ldi[(X-x/2)/Ai]) g , 

where again d = a, b. As exp ected , only electrons with opposite spin tend to form pairs. 
Inserting Eq. 4.35 into Eq. 4.39| , we obtain three averages. In the relative charge sector, 



(4.39) 



ii^exp y ±{ — 
±«4~(A X ), 



exp ( ± i^- (x/2) + Kb {-X/2) + 9^ (x/2) 9%( X /2)j 



Kb "(X/2) + Kb (-X/2) + 89"- (x/2) - 86% (x/2)] 



(4.40) 



where \ — X /A£, and we have without loss of generality chosen X — 0, since the average is independent of X by 
translational invariance. The function g%(x) satisfies 



< <&- (0) < 1, g% ~ C ie - C « X » 1, 



(4.41) 



with C\ and C2 order one constants. The exponential decay arises from the separation of the two 6% fields, whose 
rapid fluctuations exponentially suppress the average. More formally, the correlator involves a "string" connecting 
x = ±x/2, which carries an action per unit length (string tension) of l/£. 
Similar reasoning leads to the results in the spin sector: 



exp ( ± i x l - [cf> d(7 (x/2) - <M-x/2) + Mx/2) + 9 da ( X /2)] 



9da(^X), 



(4.42) 



where gda(x) are functions with the same properties as g ab p , Eq. 4.41 



The fourth factor emerging from the averages in 



Eq. 4.39 is the exponential of the 4> ab field, which is not 



averaged over. At this point, the refore, th e (^ d ±}g fields 
are still operators. Using Eqs. 4.40-4.42, they may be 
cast into the form 

{% d± (X,x)) g = ±A d (x) e^ x \i mmi ), (4.43) 



x % £ (enforced by g~ b p and g da ) to neglect the x de- 
pendence of <f>. The overall sign arises from re-ordering 
the Klein factors. Physically, we may now interpret <j) 
as the usual U(l) phase of the superconducting order 
parameter. The prefactor is what is conventionally 
interpreted as the pair wavefuncion in a superconductor, 
and has the form 



where (f>(X) = \/^4> ab {X/ A£) , and we have used the fact 
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A aO)= 2^9ab i X /09a*(x/0, 
Z7TJ; 



(4.44) 
(4.45) 



The relative minus sign between (^+)gi (&d_) g implies 
that the pair wavefunction is a spin singlet, as is demon- 
strated by rewriting this result as 

$maLjp(X,x)) g = 8 lJ A i {x){5 a] 8 Pi - 5 ai 5^)e i<t>{x) . 

(4.46) 

As in conventional superconductors, the ground state of 
the system is a singlet and, therefore, SU(2) invariant. 

The relative sign between ($ a+ ) g , ('&b+)g indicates 
that the pair wavefunction has d-wave symmetry in mo- 
mentum space. This is illustrated in Figs. [6-10]. The 
precise nature of the wavefunction, i.e the distinction be- 
tween d x 2_ y 2 and d xy pairing, depends on the positions 
of (a, b) on the Fermi surface. To emphasize this point, 
we now calculate the pair wavefunction in real space. 

The most general pairing operator in coordinate space 

is, 



*(R, r, a, 0) = ip a (R + ~)^(R - |) (4.47) 

where R = (X, Y) is the coordinate of the center of mass, 
and r = (x, y) is the relative distance between the pairing 
electrons. Here a, /3 are the spin indices of the electron 
pair. 

The pair fields * (Eq. |4.32j ) and * (Eq. fi47b are es- 
sentially related by a Fourier transform in the transverse 
(P, i) indices. To make this explicit, we must keep track 
of boundary conditions. For OBCs, the right and left 
movers are standing waves in the transverse direction, 
and 

xsm(k yl y), (4.48) 

where k V i are the transverse momenta defined in Eq. [2.8| . 
For PBCs, because the system is translational invariant, 
the decomposition is the usual Fourier one, 



V>a(r) ~ Vffio(x)e lk - r + i^ Lla {x)e-^\ 
where the fermi vector kp^ = (Jepi,2iri/N). 



(4.49) 



Consider first OBCs. Using Eqs. |4.47f44.48| , 

(*(R,r, a,/3)) g ~ ^ I sia(k yi yi) sin(k yj y2)tpi a {X + x/2)ipjp(X - x/2) 

^Y^sm{k yiyi ) S m{k yj y 2 ) RiaLjf3 (X,x)) e *rtx+*/*)-*F t &-'/*) 
ij L 9 



(4.50) 



In the second line, we have used the fact that the non-vanishing ^ operators pair right and left moving fermions, and 
have thereby dropped ^rr and ^ ll contributions. Using ^ LiaRjpiX, x) = —^ R j0Lia(X,—x), and Aj(— x) = Ai(x) 
and Eq. 4.46, Eq. 4.5C leads to the familiar form 



(*(R,r,a,/3)) g = $ d (y,r) Xa/3 e^. 
The spatial and spin parts of the Cooper pair wave function are, 

$£ BC (F,r)= 2A i (^)cos(fc i r i a;)sin(fc 9i y 1 )sin(fc 2/i y 2 ), 

i— a, b 

Xaf3 = Sa^Sfji - SalSpt- 



(4.51) 



(4.52) 



The positions of the electrons are denoted as 3/1,2 = (Y ± |). Because of the hard-wall boundary conditions, the 
spatial part of wave function depends on the transverse center-of-mass coordinate Y. Because A a A{, < 0, $° (r) 
has d-wave symmetry in real space. 



For PBCs, the results are quite similar. The pairing 
operator retains the same form of Eq. 1.51, with instead 



tive coordinate r because the system is translational in- 
variant. Once again, the symmetry is d-wave-like. 



*rf BC (r)= Y, 2Ai(:r)cos(k Fi r). 

i— a,fe 



(4.53) 

In this case, the wave function only depends on the rela- 
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D. Even chain PBCs 

Turning to the case of even N with PBCs, the generic 
presence of transverse umklapp (u) interactions allows 
for a new RG instability. This possibility is realized 
for N = 4, as we have found by numerically integrat- 
ing the extended equations of appendix B. In certain 
regions of the phase diagram, the Cooper (c) couplings 
become asymptotically irrelevant, and instead the trans- 
verse umklapp and forward scattering interactions be- 
come dominant. As before, this occurs in only two bands, 
which we will again denote a and b. In order for these 
two bands to be connected by two-particle umklapp pro- 
cesses, they must satisfy \a — b\ = N/2. More careful 
attention to the numerics shows that, in addition, these 
relevant couplings satisfy 



7 u i&>°> 



/a"a<0,A CT S <0 



(4.54) 
(4.55) 



Eqs. |4.54| - fl55| appear quite similar to Eqs. [4.20| ~ |4.2l| , 

already encountered in the generic case. In fact, care- 
ful study shows that the instability encountered here 
is mathematically equivalent, after a relabeling of the 
bands, to the earlier case. Instead of repeating the anal- 
ysis of the previous subsection ad nauseum, we will there- 
fore instead only sketch the essential points of the parallel 
treatment needed here. 

To account for the change in paired bands, we combine 
the chiral boson modes into the modified canonically con- 
jugate fields 



(4.56) 



Defining spin and charge bosons as in Eq. 4.14, the in- 
teraction terms become 



(i) 



i— a : b 

if „„„r. rrzip- 



n 



+4u L a l cos( V47r<^ 6 ) sm(V27r0 o(r ) sin(V27T0, 



1 ba) 



(4.57) 



where < jf nh = ± (\> p ^)j\f 7 2. This is of the same form 



as Eq. 4.19, and the semiclassical analysis is identical, 
with 6^6 and <f> <-> <fi. All the subsequent steps of the 
analysis carry through with small modifications. The 
non-vanishing partially-averaged pair-fields expressed in 
terms of band indices are 



(4.58) 



where <f>(X) = ^pK<\> ah (X/A£), and the gap functions in 
momentum space are 



a 6 = -A^-(z/^ 6ct (:c /£). 



(4.59) 



Notice that electrons with opposite k x but the same k y 
are paired. This implies that the Cooper pair carries 
nonzero transverse quasi-momentum. To clarify the situ- 
ation further, we now specialize to the case of primary in- 
terest, N = 4. Following the previous subsection, w e can 
again put the pairing operator into real space. Eq. 4.51 
continues to hold, but with 



®N=i(Y,x) = ^2 2A. i (a;)cos(fc Fi a;)( 



-4|A„|(x) cos(k Fa x) sin(7rF), (4.60) 



where we have u sed r eflection symmetry which implies 
g aa = g ba . in Eq. 4.59 , and hence A a = — At, as well as 
kp a = kpb = tt/2 in this case. Note that Y takes on 
integer and half-integer values, so that <i> is real but can 
vary in sign. If one imagines wrapping the four chains 
around into a cylinder, the Y dependence is simply a su- 
perposition of the m = ±2 angular momentum states, 



sin(Trr) = sin(20). 



(4.61) 



where Q = ttY/2 is the angle around the cylinder. For 
this reason, we call this a Cylindrically Extended (CEX) 
d-wave state. Note that since the superposition here is 
purely real, the state does not carry any spontaneous 
current. 



V. PHASE DIAGRAMS 



In the previous sections, we have described the RG and 
bosonization technology necessary to analyze a weakly- 
interacting one-dimensional Fermi system for any generic 
set of parameter values. We have, of course, applied these 
methods to study the particularly interesting case of the 
iV-chain Hubbard models. The detailed calculations in- 
volve lengthy but straightforward numerical integrations 
of the RG equations and mapping out the ensuing pair- 
ing instabilities as a function of N, n, and t±/t. The pri- 
mary results of this work are the phase diagrams shown 
in Figs. [6-10]. For the most part these stand on their 
own, but we will comment on a few points. 
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FIG. 6: Phase diagram of the two-chain Hubbard model 
with weak repulsive interactions. 



FIG. 8: Phase diagram of the 3-chain Hubbard model for 
PBCs. 
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FIG. 7: Phase diagram of the 3-chain Hubbard model for FIG. 9: Phase diagram of the 4-chain Hubbard model for 
OBCs. OBCs. 
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FIG. 10: Phase diagram of the 4-chain Hubbard model 
for PBCs. The states marked with an asterisk exhibit 
CEX pairing. 



A. Commonalities 

1. Band transitions 

In the weak-coupling limit, it is natural that the gross 
features of the phase diagrams are dictated by the non- 
interacting band structure. In particular, the n—t±/t 
plane is divided into distinct regions, in each of which 
a particular number Nf of Id bands are partially filled 
(and hence not inert). The boundaries between these re- 
gions constitute band transitions, which generally survive 
as phase boundaries in the weakly interacting system. 

While the band transitions need not be the only phase 
boundaries in the interacting system, they usually form 
the most noticable divisions of t he phase space. To lo- 
cate them, one must solve Eqs. ( |2.7| |2.8| , 2.10- ^.11 ) for 
the lines along which e a (0) = /i (band a is just empty) or 
£a(7r) = fJ- (band a is just full), for each a. These curves 
are indicated in the figures by heavy lines. The shad- 
ing of the regions separated by the heavy lines indicates 
the number Nf of partially filled bands, with white corre- 
sponding to Nf = 1 and the darkest shade corresponding 



most empty or almost filled. This gives rise to a very 
small Fermi velocity, v E C t, in the weak-coupling RG. 
For such a small velocity, the dimensionless couplings 
acting purely within the band E are greatly enhanced: 
See = 9ee/{^ve) > 9aE, for a ^ E. Physically, this 
increased scattering is due simply to the large density of 
states near the Id van Hove singularity at the edge of the 
band. 

Strictly speaking, the RG equations as we have derived 
them are not valid directly at the band transition. This 
is because the spectrum of an empty /full band is not 
relativistic but quadratic (u> ~ k 2 /2m), demanding a dif- 
ferent (anisotropic) scaling. We can, however, approach 
very close to the band transition in the weak-coupling 
limit. That is, provided that we keep U/ve *C 1 (a much 
more stringent requirement than U/t <C 1), Eqs. 3.7-3.1C 
remain valid. 

In this region, the nearly vanishing Fermi velocity pro- 
vides a useful small parameter. Indeed, those couplings 
within band E are (at least initially) much larger than 
those involving any of the non-critical bands. For the 
initial stages of integration of the RG flows, then, these 
"edge-band" interactions dominate the evolution of the 
couplings. In particular, they lead to a decoupling of the 
critical band from the remaining degrees of freedom. To 
see this, consider the evolution of the interactions be- 
tween band E and another band i. In the Cooper chan- 
nel, 



-IE 



'EE 

"-EE <-EE 



W C EE 



2 C EE 



Because the initial value for the Hubbard model, 



iE 



(0) 



<&(0), 



(5.1) 



(5.2) 



hap pen s to be one of the eigenvectors of the matrix in 
Eq. 5.1, the solution is particularly simple: 



cf B (0 = < E (0 * Cie 



.£2.1 



(5.3) 



where C%, Ci are constants of 0{U). 

This exponential decay holds provided the running 
couplings within band E {cf EE , c EE ) remain large com- 
pared to the inter-band couplings. Examination of 
Eqs. 3.9-3.10 shows that the intra-edge-band couplings 
relax logarithmically to zero in the initial flow regime. 
The above equations are thus valid for 



to N f 



N. 



pa p 
C EE' C EE ~ T c kEi< 



■kE 



U 

Vk 



(5.4) 



2. Band edge phases 

On one side of such a band transition, the "criti- 
cal" band (let us denote its index by a — E) is al- 



Following Eq. [T| out to this, we see that the couplings 
to the critical band are indeed exponentially suppressed 



cfe(i = v k /U) = 4cfe(l = v k /U) 



C\ exp[-C 2 vk/vE], 
(5.5) 
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where C\ and Ci are order one constants. 

To complete the argument that band E becomes de- 
coupled from the others, we must now show that this 
suppression persists into the regime of divergence of other 
couplings, i.e. to the cut-off scale I*. To do so, we employ 
the alge braic relations valid in the asymptotic regime 

(Eq. |4.8| ). Consider, for example, the relation derived 

from Eq. |3.9| . The contributions on the right hand side 
can be separated into singular and non-singular parts, 



-OLki,EG^ E Ci E + (non-singular terms), (5.6) 



where k,l ^ E. Here we have used the relaton Cfj = 
as usual. Since the factor ctki.E ~ 1/Ve is singular when 
ve — > 0, there are only two possible options. Either the 
singularity in au,E is cancelled, and Cue ~ y/vE, or the 
singular and non-singular parts equal zero separately and 



GIe = o. 



(5.7) 



The former possibility i s inc onsistent with the exponen- 
tial sup pression in Eq. 5.5, so we expect that instead 
Eq. p/ij holds and the couplings c^e will flow to zero in 
the asymptotic regime. This is indeed observed in all nu- 
merical integrations of the full RG equations near a band 
transition. 

This decoupling implies that the low energy structure 
of the system is obtained by adding the single gapless 
charge and spin modes of the critical band to the low 
energy structure of the remaining bands that would have 
occurred were band E inert. To determine the phase of 
the Hubbard model in the band-edge regime, therefore, 
we may simply add C1S1 to the gapless mode content 
on the other side of the band transition, in which the 
critical band is indeed inert. If the phase on this side is 
CnSm, then the band-edge result on the other side of 
the transition line is 



CnSm ~\~ C1S1 — Cn+iS 



+iSm+i- (5. 



B. Specific features 

1. "d-wave" pairing 

Probably the most striking aspect of the phase dia- 
grams is the ubiquity of paired states - i.e. gapping out 
of the spin modes in at least some of the bands. Follow- 
ing the methods of section IV, these pairing instabilities 
can be associated with a gap function defined at the al- 
lowed discrete points on the 2d Fermi surface. Except 
in certain regions of phase space in the four-chain model, 
this gap function has an approximate "d-wave" form. For 
N > 3, a d-wave gap has an interesting consequence in 
this context: the discrete transverse wavevectors k ya can 
coincide with the nodes in the pair wavefunction. This 
indeed occurs, e.g. near half-filling for N = 3 with OBCs 



for tj_ < y/2t, giving rise to simultaneous dominant su- 
perconducting correlations and power-law antiferromag- 
netism. The d-wave interpretation begins to break down, 
however, for N — 4 with OBCs, where several gapless 
spin modes are present for small tj_/t. The unusual dis- 
tribution of gapped and ungapped modes on the Fermi 
surface in this case is, we expect, a consequence of the 
one-dimensional weak-coupling limit taken here, which 
becomes rather restrictive for larger N. A more com- 
plete discussion of the approach to a two-dimensional 
weak-coupling limit is described in section VI. 



2. Peculiarities of PBCs 

The systems with PBCs exhibit a number of (theoret- 
ically!) interesting peculiarities. For N ^-3, some of 
these have been pointed out by Arrigoni,EH who per- 
formed a similar weak-coupling analysis. For odd N , the 
effects of PBCs can be expected to be rather severe. In 
the Hubbard model, they break particle-hole symmetry, 
which has the effect of eliminating reflection symmetry of 
the phase diagram around half-filling. Furthermore, half- 
filling no longer coincides with the conditions needed for 
umklapp processes at the Fermi level. This has the happy 
consequence that our generic treatment (which ignores 
these umklapp interactions) remains valid at n = 1 . One 
then finds the rather surprising result that the half-filled 
system has gapless charge excitations and a spin gap, 
precisely the opposite of what is expected in the strong- 
coupling limit with OBCs, where there is a charge gap 
(~ U), and the effective Heisenberg model (with odd N) 
is expected to have a gapless spin mode. Some partial 
understanding can be gained from the fact that an odd- 
chain Heisenberg model with PBCs is frustrated, and can 
be shown (at least for N = 3) to indeed have a spin gap. 
However, the absence of a charge gap is a weak-coupling 
result, and indicates the existence of a metal-insulator 
transition at half-filling as U is increased. Interestingly, 
the weak-coupling paired state can be either of d xy or 
d x 2_ y 2 type, as indicated in Fig. [8]. 

A different sort of feature arises for PBCs with N = 4. 
Although this situation retains particle-hole symmetry, 
there nevertheless exist regions (the largest occurs for 
weak doping with t±_ > t) in which the finite transverse 
size has a severe effect. This CEX d-wave phase has a pair 
wavefunction whose phase depends upon the transverse 
center-of-mass coordinate of the pair. This phenomena is 
described in section IVD, and is certainly special to the 
one-dimensional cylindrical geometry considered here. 



3. Extreme asymptotic instability of the C1S0 phase 

A number of authors have predicted the existence of 
a CI SO paired state for a weakly interacting two-chain 
Hubbard ladder. As noted (in proof) in Ref. O in the 
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asymptotic limit U/t — > + , the CI SO phase in fact oc- 
curs only for infinitesimal doping, being replaced every- 
where else by the C2S1 phase. For reasonable (but still 
small) values in the range 10~ 6 < U/t < 1, however, the 
C1S0 phase still appears as dominant. 

To understand this result requires a more detailed ex- 
amination of the asymptotic regime near the RG diver- 
gence. To do so, we again consider the algebraic relations 
described in Sec. 



IV A 



Using Eq. 4.£ , we look for a solution of the resulting al- 
gebraic equations for which Cf 2 7^ 0, corresponding to a 
C1S0 phase. Some straightforward calculations give the 
unique answer 

1 



C?i = C° 22 = --(1 + y/l - 4au, 2 Cf 2 ), (5.9) 
-(1 + ^,2X^=0. (5.10) 
Now consider the stability of this solution. From 



Eq. 3.10, the difference between c\ x and c 22 obeys 



■Ki-42)Wi + 4 2 )- (5-ii) 



Since are negative near the divergent point I*, the 
difference of Cooper couplings d{\ — c 22 is relevant, and 
the above solution is in fact unstable. For the special 
initial value c\ x — c 22 = 0, which is attained in the limit 
n — ► 1, the system is specially tuned to an unstable equi- 
librium, and the C1S0 phase enjoys a small region of 
existence. Moving away from half-filling, however, the 
inequality of Fermi velocities in bands 1 and 2 destroys 
this fine-tuning, driving the system away from the C1S0 
state. In the asymptotic limit, the stable solution is in 
fact the much simpler C25*l flow. This point has been 
missed in other calculations, owing to the assumption 
of equal Fermi velocitiesEjO and the lack of « jCareful 
stability analysis of the asymptotic regimes. Elf EH We 
emphasize, however, that numerically this instability is 
extremely weak. For even relatively weak couplings with 
10 -6 < U/t < 1, we find that the C1S0 phase remains 
quasi-stable, occupying in fact the majority of the two- 
chain phase diagram. 



VI. DIMENSIONAL CROSSOVER 

In this section, we discuss how the system approaches 
2d behavior as N — ► 00. The limit is actually quite sub- 
tle, and we will consider two distinct ways of performing 
it. The simplest procedure is simply to attempt to pre - 
serve the validity of the RG as presented here (Eqs. 3.7- 
3.10). This Id weak coupling limit can be realized in 
principle for any fixed (but large) N for sufficiently small 
U , but taking N — > 00 actually requires that the in- 
teractions vanish as well. A more physically appealing 
approach is the true 2d weak coupling limit, in which 
N — > 00 for a fixed (but small) U. In this case, the RG 



as constructed so far must in principle be supplemented 
by additional interactions. 



A. One-dimensional weak coupling limit 

We fir s t co nsider the naive limit of the RG flows 
(Eqs. 3.7-3TC) as N — > 00. Recalling the results of sec- 
tion II (Eq. |2.46|), to retain the validity of these equa- 



tions, the interactions must be simultaneously taken to 
zero, with U ^ t/\nN. Since this constraint is only log- 
arithmic in N, this is actually not a strong restriction 
even for reasonably large values of N, and might indeed 
be physically relevant in some systems. 

In the large N limit, the RG flows are dominated by 
those terms involving sums over intermediate band in- 
dices; which effectively increase these terms by a factor 
of AT. To make the largest terms of order one for large 
N, we introduce the rescaled coupling constants 



p = 1 



(6.1) 



and similarly for the other interaction channels. Note 

that, with Hubbard initial conditions, the values of cf,- 

— j 



arc o f an order of one (in TV). Inserting this into Eqs. 3/7 
3.10 and dropping the 0(1 /N) terms, one finds that 
the forward and umklapp scattering vertices are exactly 
marginal (unrenormalized) . The Cooper channel inter- 
actions obey the simplified equations 



k 



(6.2) 



= -Jj J2^k £ k 3 + clkKj + 7}cf k c%j). (6.3) 



The Hubbard model initial condition that 4cf-(0) 



cfJO) is preserved by Eq. |6.3| , so they may be collapsed 



into the single non-trivial RG equation 



(U: 



dl 



k 



(6.4) 



Note that the right-hand-side of Eq. 3.4 has the form 
of matrix multiplication. This implies that, provided the 
initial couplings matrix is diagonalizable, each eigenvalue 
Xi of the matrix d 7 evolves independently according to 



d\j 
~~dT 



1 

'N' 



(6.5) 



Focusing for simplicity on PBCs, the initial value of c° 
is 



^■(o) = (-) 



u , 1 



2n y<ViVj 



(6.6) 



Since this has the form of an outer product, it is propor- 
tional to a projection operator onto the vector 1/ ^/v[. 
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It thus has Nf — 1 null vectors ancLa-single non-trivial 
eigenvector proportional to l/y^.OCJ The eigenvalues 
are 



A,-(0) 



v /TT ' ' 



' 2ti ^ v k 



(6.7) 



From Eq. 6.5, the Nf zero eigenvalues are unchanged 
under the RG, while Ai obeys 



MO 



ATAi(O) 
N + Ai(0)f 



(6.8) 



where A x (0) = [U/2n) £ fc l / v k- Since Ai(0) > 0, it is 
marginally irrelevant and flows to zero. This implies that 
all the Cooper interactions flow logarithmically to zero. 

To connect with previous two-dimensional treatments, 
we may define a simple continuum limit: 



N ^ 



V(k 
77 dk 



h 



i), 



2tt 



(6.9) 
(6.10) 



which gives the RG equation 



as 



= - y ^n*,,^^.^)- (6-ii) 

derived previously by ShankarEl directly in 2d. 



vertex. Next, note that, in weak coupling, the range of 
momentum shifts is highly constrained: 



\Ak v \ = 2ir\8\/N < 2TT6 max /N ~ At ± /t, 



|AA)y| ^ |Afcy| max 



tx - 
= TT — e 

t 



const. £/£/ 



(6.12) 



as can be seen from Eqs. |2.44| and |2.15| . For any non- 
singular interaction, the initial coupling constants are 
reasonably smooth in momentum space, so that all the 
shifted interactions in the narrow range \k y s\ < |Afe y | are 
essentially equal in magnitude. Since each shifted vertex 
then has the same initial conditions and obeys the same 
RG equation as an unshiftcd coupling, it remains so un- 
der the RG, and the original equations remain sufficient 
to study their evolution. 
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FIG. 11: An example of renormalization of from 
shifted vertices in Eq. 6.13. 



B. Two-dimensional weak-coupling limit 

On reflection, the agreement with approaches directly 
in two dimensions is perhaps surprising, since the RG 
equations used above are valid only for U/t < 1/lniV 
as JV — ► Co. To study the true 2d weak-coupling limit 
(with U/t -C 1 but fixed as N — > oo) requires con- 
sideration of the additional shifted interactions (such 
as 7i s ) introduced in section IIC. Fortunately, one can 
show that, even upon including these interactions, the 
modifications of the RG equations are actually negligible 
in weak-coupling. Rather than belabor this reasoning, 
which is essentially discussed already in, e.g. Shankar's 
review articleJ23 we will only schematically indicate how 
this comes about. 

Once the additional shifted interactions are included 
in the RG, we must worry about two questions. How 
do these new vertices renormalize, and how do they feed 
back into the flow equations for the unshiftcd couplings? 
In answer to the first question, under normal conditions, 
the shifted interactions renormalize almost identically to 
their unshifted counterparts, at least in the initial stages 
of the RG. This is because for each process involving 
two unshifted vertices feeding into an unshifted vertex, 
there is an analogous process involving the same vertices 
shifted, usually feeding back into the analogous shifted 



It remains to answer the second question. There are 
a few additional processes involving the shifted interac- 
tions, such as the one shown in Fig. [11], which feed back 
into the original Cooper and forward scattering channels. 
The feedback into the forward scattering channel is neg- 
ligible for the same phase-space reasons that render them 
exactly marginal for N — > oo above. More subtle is the 
feedback into the unshifted Cooper channel. Once shifted 
vertices are included, an intermediate sum over 5 allows 
for non- vanishing contributions of the form 

^ = -" + ^£W«)*--) (*)+•". (6-13) 

8 

in which the (shifted) forward scattering vertices feed 
back into the Cooper channel even as N — > oo (since 
8max./N is finite in this limit). However, the phase 
space for these renormalizations is considerably smaller 
than the processes already included. Roughly speak- 
ing, these additional terms are down by a factor of 
A/ir r~j exp(— const. t/U) from the others, since the range 
of angles of the intermediate momenta are restricted to a 
width |Afc y | max . <C 2tt. Furthermore, the number of al- 
lowed terms continues to decrease as the RG proceeds to 
lower energies, since the band curvature t±/t effectively 
grows under rescaling. 
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C. Instabilities for 1 < N < oo 

In the large N limit, therefore, additional (shifted) in- 
teractions are present at weak but finite couplings, but 
do not modify the RG equations or their analysis as 
presented above. ,-^t TV = oo, then, the 2d metal is 
marginally stableO What occurs in weak coupling for 
large but finite N? Our numerical results suggest that 
instabilities always persist when feedback of forward- 
scattering interactions into the Cooper channel is in- 
cluded. However, this cross-coupling is an 0(1/N) ef- 
fect, and sois extremely weak for large N. Specifically, 
these terms can only begin to affect the flows once the 
Cooper interactions have themselves renormalized down 
to order 1/N. Since the N — oo flows are logarithmic, 
this occurs only after a rescaling b = e l ~ e N , so that the 
characteristic energy gaps (and critical temperatures) of 
any paired states should obey 

A N <e~ N , (6.14) 

with a prefactor which is not determinable by such coarse 
arguments. 

This exponential decrease of the energy scale for pair- 
ing is a signature of the rather robust stability of the 
generic Fermi liquid. This weak-coupling result, however, 
does not make any statement about pairing instabilities 
for strong repulsive interactions. Non-generic situations 
can, of course, give rise to much larger energy scales, even 
at weak coupling. Of particular importance in highly 
anisotropic repulsively interacting systems is the spin- 
density-wave instability. Because this requires nesting in 
weak coupling, the associated interaction vertices have 
been thrown out in our calculations. 



VII. DISCUSSION 

The principal results of this paper are the iV-chain 
weak coupling phase diagrams, described in detail for 
small and large N in sections V and VI, respectively. 
We now conclude with a discussion of the implications of 
these results for both ideal finite N systems (accessible 
via numerical calculations) and for true quasi-ld systems, 
where non-zero inter-ladder couplings need to be taken 
into account. 



A. Numerics 

Numerical calculations have the advantage that direct 
comparisons with Id models can be made. Recall that 
each phase is characterized at the simplest level by its 
number of gapless charge and spin modes. These num- 
bers can be measured numerically in a number of ways. 
Most directly, the lowest lying charge and spin excitation 
energies can be determining by comparing ground state 



energies (in, e.g. exact diagonalization or density-matrix 
RG methods) with particles added or spins flipped. Such 
measurements can also be refined to determine the ener- 
gies for the lowest-lying excitations with a definite parity, 
which can be related to the band indices used here. The 
total number of gapless modes (both charge and spin) can 
in principle be extracted alternatively from the coefficient 
of 1 JL (where L is the chain length) in the finite-size cor- 
rection to the ground-state energy density. 

The parity of the ground state and low-lying excited 
states are also accessible in weak coupling. Here we focus 
on the C1S0 phase in a 2-chain system as an example. 
The parity operator of a 2-chain system is 

P = exp(t7rJV 2 ), (7.1) 

where N2 is the total number of particles in the anti- 
bonding band. The parity operator P commutes with 
the Hamiltonian, so that parity is a good quantum num- 
ber. If the ground state is a linear superposition of states 
with odd/even N2, it has odd/even parity. In the CI SO 
phase, if total number of particles is even, N 2 is even 
because electrons pair up in both bands. Therefore, the 
ground state has even parity. However, if the total num- 
ber is odd, the analysis is complicated and remains an 
open question for further study. The parities of excited 
states is determined by commutation relation between 
the corresponding creation operators and parity opera- 
tor. The Bosonic field operator <^"(p), which is the only 
gapless mode in the C1S0 phase, creates a density exci- 
tation with momenta p. If we express N2 in terms of the 
Bosonic fields 

N 2 = VMM 00 ) - M-oo)}, (7.2) 

it is easy to show that <p p ^(p) commutes with the par- 
ity operator. This implies that the excited state has the 
same parity as the ground state. In other words, the 
excitation carries even parity. Consequently, numerical 
calculations should find a charge gap in the odd parity 
channel, despite the exsistence of a gapless charge mode 
with even parity. A simple test is the correlation function 

C{x) = (Ap(z)Ap(O)), (7.3) 

where Ap(x) = c\ a c\ a (x) — c\ a C2 a (x). Since P^ ApP = 
—A/?, C(x) has only odd-parity (relative to the ground 
state) intermediate states, and should therefore decay ex- 
ponentially, C{x) ~ e - A p- x / v \ 

It is often more convenient to compute correlation 
functions rather than ground state energies (indeed, in 
a Monte Carlo calculation, this is essentially the only op- 
tion). In this case, the information that can be most 
reliably assayed is the presence or absence of charge and 
spin gaps. If there is a true gap in either sector, the corre- 
sponding correlators are expected to decay exponentially 
in space. To probe the charge sector, the correlators of 
interest are those of the density p — c^ a c a and the pair 
field A = CfCj_. In the spin sector, the corresponding 
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operator is simply the spin S = cJ^Cg. In principle, a 
detailed examination of the Fourier content of the cor- 
relations should identify the Fermi momenta of the gap- 
less spin and charge modes (even more information than 
their number), but this is quite difficult in practice due 
to finite-size limitations imposed by the numerics. 



t',U' 



1. Weak interactions 



To proceed, let us imagine repeating the weak-coupling 
RG with these additional interactions. Upon first inte- 
grating out the large wavevector modes ( | fc^ | > A) on 
each ladder, t' and U' will suffer small renormalizations, 
and other new interactions will also be generated. Two 
of these are of particular importance: two-particle hop- 
ping processes, in which two Fermions are simultaneously 
transferred from one ladder to a neighbor, and four- 
particle "pair-density" interactions, in which Fermions 
interact energetically on neighboring ladders in a man- 
ner quartic in the density, but no charge is transferred . 
These will occur with a pair-hopping amplitude t" and a 
pair-density interaction U" , which are approximately 



t" 
U" 



(tru 2 /t\ 

(U') 2 U 2 /t\ 



(7.4) 
(7.5) 



FIG. 12: A 2d array of Hubbard ladders, with weak inter- 
ladder hopping t' and inter-ladder density-density inter- 
action U' . 



B. Experimental consequences 



Comparison with experiments is more challenging. In 
particular, it is inevitably the case that in any candidate 
compound there is at least some residual coupling be- 
tween ladders as shown in Fig. 12. In this sense, all real 
materials are at best quasi-ld. It is to the features of 
such quasi-ld ladder materials and the regime of validity 
of the previous results in this context to which we now 
turn. The discussion will be kept at a general level, using 
only scaling considerations. We take as a model a two- or 
three-dimensional regular array of A^-chain ladders, the 
precise geometry of which is not crucial, although spe- 
cial cases resulting in Fermi-surface nesting will not be 
addressed. A microscopic electronic model for such an 
array would involve the hopping amplitudes and interac- 
tions both on and between the ladders. We will assume 
that the former are of the Hubbard form studied in the 
previous sections. The latter generically introduce two 
new energy scales: an inter-ladder hopping amplitude t! 
and an inter-ladder density-density interaction U'. For 
now, we assume that, at least t'/t, U'/t <C 1. 



in the weak-coupling limit (see Fig. [13]). For general- 
ity, we shall keep t" and U" as independent parameters. 
Other interactions are of course also generated, but are 
either of similar type but much smaller magnitude than 
those already considered, or are higher-order and hence 
at least perturbatively irrelevant. 




(a) 



(b) 



FIG. 13: Lowest-order diagrams illustrating the genera- 
tion of inter-chain pair-hopping t" and pair-density in- 
teractions U" from the bare inter-chain hopping t' and 
density interaction U'. 

At this point, we proceed with the RG as before, work- 
ing perturbatively in t'/t, t" /t, U'/t, and U"/t. Like 
the original ladder parameters, these will also rescale 
and non-linearly renormalize t hemselve s and other cou- 
plings. The corrections to Eqs. 3.7- 3.10 will, however, be 
negligible provided the values of the running couplings 
t'/t\ h t"/t\i, U'/t\u U"/t\i < g t (l). Initially, this of course 
requires t' ,t" ,U' ,U" -C U, but the constraints become 
stronger as we iterate to lower energy scales. In partic- 
ular, for the divergences encountered purely within the 
ladder RG to be essentially unchanged requires that the 
running inter-ladder couplings be negligible compared to 
one (since the relevant ladder couplings become of order 
one) at the scale I* ~ ct/U , where c is a constant. Simple 
power counting gives 
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(7.6) 
(7.7) 
(7.8) 
(7.9) 



since t' , t" , U' and U" represent two, four, four, and 
eight- Fermion operators. The hopping if /t thus grows 
exponentially, the two-particle processes t" /t,U' jt scale 
only logarithmically (are marginal), and the pair-density 
interaction U" is strongly irrelevant. The Id RG results 
are then valid up to the instability scale I* provided 

t' < te~ ct / u ~ A, (7.10) 

t" < U, (7.11) 

U' < U, (7.12) 

U" < te 2ct ' u ~ t 3 /A 2 , (7.13) 

where A is the energy scale of the gap in the Id system 
(note that factors of (t/U) in the prefactor are not cap- 



tured within a one-loop RG treatment, so this is a rough 
estimate and not a strict asymptotic statement). The re- 
quirement of t' < A has a simple physical interpretation: 
for tf > A, it is favorable for singlet pairs to break up to 
reduce their kineti c energy, d estroying the paired state. 
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Provided Eqs. |7.10| - |7.13| are satisfied, the strong- 
coupling analysis of section IV holds, and the paired 
bands are adequately described by the single collective 
phase mode <fi p+ (and its conjugate 9 P+ ). For concrete- 
ness, we now specialize to the two-chain case, where there 
are no additional bands. The single-particle tunneling 
operator conjugate to t' then involves exponentials of the 
dual fields 4> aa (c.f. Eq. 4.28), which fluctuate wildly and 
are exponentially suppressed (strongly irrelevant). Simi- 
larly, the density-density interaction U' is also negligible 
due to strong fluctuations of the relative-displacement 
mode 9 P ~ . The remaining two couplings (t" and U") sur- 
vive, and have simple interpretations. The pair-hopping, 
t" simply hops a single Boson between neighboring lad- 
ders, and is hence like a Josephson coupling. The pair- 
density coupling U" is effectively a density-density inter- 
action between Bosons on neighboring ladders, which are 
created by the pair fields of section IV. In terms of the 
phases, the effective Hamiltonian is 



E 

(nn') 



t" cos V^(r n + + ^ cos 2v^(0£ + - O 



(7.14) 



where the index n labels the different ladders, and we 
have taken the cut-off scale A = 1 for simplicity. The fac- 
tor of £ -2 in the pair-density interaction reflects its irrele- 
vance at the non-interacting Fermion fix ed point, and re- 
sults directly from integration of Eq. 7.9. As discussed in 
Rcf. |l3| the stiffness K is not exactl determinable within 
the weak-coupling RG for generic parameters. However, 
various arguments suggest K > 1/2^3 and in particular 
K w 1 close to half-filling in the two-leg (N = 2) ladderfed 
at weak coupling. 

To determine the nature of the ground state of the sys- 
tem, we must now address the physics at energy scales 
below A. To do this, we use a different RG rescaling, 
using standard sine-Gordon techniques e3 This gives the 
linear flow equations 



^(s-ijf+or) 2 ,^ 2 ], 

d 



(7.15) 



jjjU" m (2 - 2K)U" + 0[{t"f, {U"f\. (7.16) 

Note that the Josephson coupling t" is relevant for K > 
1/4, while the pair-density interactions are relevant for 
K < 1. There is therefore no region of stability for 
truly one-dimensional behavior, regardless of K. Most 
probably 1/2 < K < 1, and both perturbations are rel- 
evant. If both the dimcnsionlcss bare interactions are 



weak, t" /t,U" /(t(, 2 ) -C 1, then the nature of the insta- 
bility is determined by the interaction which renormal- 
izes to large (order one) values first. Simple algebra thus 
predicts that pair-tunneling dominates for 



f „^2-2K fjjn\ 2-1/(2*0 



(7.17) 



and the ladders phase lock into a bulk superconducting 
(SC) state. In the opposite limit, pair-density interac- 
tions dominate and lead to a paired-insulator or CDW 
state. Note that for K close to 1, the SC state domi- 
nates for all but extremely small t" . Ge neric ally, though, 
as t" is reduced below the limit of Eq. 7.17, the system 
makes a transition to the CDW state. 



Let us now use the estimates in Eqs. 7.4-7.5 to deter- 



mine the bulk phase diagram for the two-chain system 
in the weak-coupling limit. As t' is decrease d, th e pair- 
ing instability occurs first, according to Eq. 7.10 , when 
t' < A. Just below this scale, it is straightforward to 
show that this instability always leads to a SC rather than 
CD W, pr ovided U/t < 1, as supposed. This is because 
Eq. 7.17 can be rewritten, using Eqs. 7.4-7.5 and the 
scaling of the dimensionless coherence length £ ~ t/A, 



t' t 
1 ~ U 



U'UA 



t 3 



(1-1/4K)/(1-K) 



(7.18) 
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and (1 - l/(AK))/{\ - K) > 1 for K > 1/2. Only for 
much smaller t' does this inequality cease to hold and 
the system go over into a CDW state. A schematic zero- 
temperature phase diagram for fixed small U jt with these 
features is shown in Fig. [14]. 



Fermi Liquid 



Superconductor 



CDW 











U' 



FIG. 14: Schematic zero temperature phase diagram in 
the t'-U' plane for an array of coupled ladders. 
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This classical model has only a single dimensionless cou- 
pling constant, as can be seen by rescaling x — > x/K. The 
superconducting transition must occur when this dimen- 
sionless value is order one, giving the critical temperature 



• „ x 2K/(AK-1) 

^ ' 7 



(7.20) 



valid within the SC region of the phase diagram away 
from the zero temperature quantum SC-CDW transition. 
Note that for t" <C t, there is a large temperature range 
T sc <C T <C Tp a i r over which the system has a "pseudo- 
gap" -like behavior. The difference between the exponent 
2K/(4K— 1) and 1/2 (the classical result for weakly cou- 
pled Id XY chains) represents a suppression of the tran- 
sition temperature due to quantum fluctuations. 



2. Beyond weak interactions 



Supposing the system is in the SC phase at zero tem- 
perature, what is the expected phenomenology as the 
temperature is varied? In weak coupling, we expect sev- 
eral large crossover ranges. For T 3> T pa i r ~ A, the sys- 
tem acts approximately in a non-interacting Id fashion, 
with small logarithmic corrections which are precursors 
of the instability to be encountered for T ~ T pa i r . Be- 
low that temperature scale pairing effectively occurs, and 
measurements (e.g. magnetic susceptibility or electron 
tunneling ) probing single-particle and spin excitations 
should exhibit activated behavior. However, although 
pairs form at T pa j r , superconducting coherence sets in 
only at a lower temperature, T sc . 

To see this, we continue rescaling after reaching I*, i.e. 
to energies below the spin gap. We can now rescale fur- 
ther, as indicated above, until the rescaled temperature 
grows to the order of the energy cut-off, at which point 
the zero-temperature RG fails. Having already rescaled 
to I* ~ ln(t/A), the effective temperature has already 
been increased to (t/A)T, but for T < A, this is still 
small relative to the cut-off t, and we can rescale further 
by the factor b ~ A/T. At this point the temperature 
is on the order of the cut-off, and thermal smearing is 
sufficient to remove any remaining quantum coherence 
at lower energies. The corresponding length L qc ~ v/T 
represents a quantum to classical crossover scale at tem- 
perature T. Fluctuations of larger size behave essentially 
classically (as can be shown explicitly by restricting the 
rescaled effective action to the zero Matsubara frequency 
modes) and can be studied using the rescaled classical 
(phase only) model 



Although the analysis of this paper has assumed U <C 
t, many two-chain models have now been convincingly 
demonstrated, numerically and in some limits analyti- 
cally, to display spin gaps when weakly doped even for 
relatively strong interactions (e.g. U > t). In fact, most 
aspects of the phenomenology in this conclusion are ex- 
pected to continue to hold even in this limit, provided 
that the inter-chain couplings are small, U',t' <C t. In 
particular, on physical grounds, we expect a transition 
to a Fermi liquid (or at least two-dimensional behavior) 
for i' > A. For t' < A, inter-chain couplings essen- 
tially never break pairs , and the physics will still be well- 
described by Eq. 7.14 at these and lower energy scales. 
Of course, in strong coupling the parameters t" and U" 
cannot be estimated usin g the w eak-coupling diagrams 
of Fig. [13] leading to Eqs. 7^-7^. In addition, it is diffi- 
cult in strong coupling to determine K: recent numerical 
simulations with U/t = 8 suggest values of K w 1/2, 
somewhat smaller (and hence less superconducting) than 
in weak-coupling. With these caveats, the remaining phe- 
nomenology should continue to hold, both for the zero- 
temperature phase diagram and for the crossovers and 
transitions at T > 0. 
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APPENDIX A: CURRENT ALGEBRA 

Current algebra methods allow, among other things, 
an algebraic calculation of the one-loop RG equations. 
Here we give a very terse description of this method. 
All currents are defined in terms of the fermion fields 
ipR/Li a (i — 1,2, ...,N), which obey the operator prod- 
ucts, 



^m<x{x,T)tjy R:j p{Q,oy 



3jj Sgf3 
2lTZi 

V^,t)^(0,0)~^ 



+ 0(1), 

+ 0(1), (Al) 



where zi = vit — ix. The operator products should be 
understood to hold when two points (x, t) and (0, 0) 
are brought close together. We therefore only need to 
keep the singular terms as replacement within correlation 
functions. As an example, consider the product JijJim- 
Performing all possible contractions gives 



Jij(x,T)J lm (0,0)~: V'L^-c " W^m/? : 

9 X. 
SimSjl + 



Am^ZiZj 



2irzi 



la 



5 



2ttz^ 
2irzj 



2ttz. 



Jli 



3 iTT 2 ZiZj 



0(1). 



(A2) 



where J denotes a 4-dimensional integral over the two 
complex planes z and w. As in any RG, we wish to in- 
tegrate out the short-scale degrees of freedom to derive 
the effective theory at long wavelengths and low energy. 
Here this is accomplished by considering the contribu- 
tions to Eq. A6 when the two points z and w are close 



together (near the cut-off scale) . We make use of the op- 
erator product expansion to integrate out the short-scale 
degrees of freedom, which gives, 



1 



TC Td 
" J R,-., J T 



Wi ) 2n{z* - w*) Rii L "- 
(A7) 



We choose a short distance cut-off a = A -1 in space, but 
none in imaginary time. For a rescaling factor 6, we must 
then perform the integral, 



1 



In 6 



lin = \ dx dr t 

la<\x\<ba J-oo {2-nYZiZ* TT{Vi + Vj) 



, (A8) 



where z is the relative coordinates in Eq. A8. The con- 
tribution to the RG equation after integration is, 



TTiVi + Vj)' 



-dl 



J Rii*J Ljj i 



(A9) 



where dl — In b is the logarithmic length scale. This 
term, when re-expone ntiat ed, renormalizes and gives 
the first term in Eq. 3.8. All other terms in the RG 



equations can be carried out by similar steps. 



We can compute the full set of operator products in a 
similar way. The ones needed are: 



Jij Jl r. 



1/2 . x , S ah Sji T 



ibc X 



4 v 27rz, 



2 v 2'KZj 



2'KZj 



5ji 



J a 



2-KZi lj ' 



2ttz. 



■Jli) 



(A3) 
(A4) 



APPENDIX B: RG EQUATIONS OF UMKLAPP 
INTERACTIONS 

For PBC, when the number of chains is even, we need 
to study the transverse umklapp interactions in addition 
to the forward and Cooper vertices. We then obtai n ad- 
ditional terms on the right hand sides of Eqs. pT? -3.1 
Denoting these by Sftj, 5cij, one finds 



where the coordinates of two operators on each left hand 
side are consecutively (x,r) and (0,0). Similar forms 
hold for the left-moving currents, but with Zi — > z*. 

The RG equations can be obtained very simply from 
the operator product expansions. We use the functional 
integral formulation which results in the Euclidean ac- 
tion, Se = J dxdrTL, and the partition function, 



^ = -%E^4(40 2 + ^) 2 } 



-{(- 2 /) 2 + ^Kf) 2 } 



+E{(^) 2 + iiKf) 2 }> 

a 



(Bl) 



Z = / 



(A5) 



To perform the RG, the exponential is expanded to 
quadratic order in Ti.. A typical term takes the form, 

\fijL I (JZiiWJijMJmWJlrnrnH), (A6) 



-{2u-jU v 



+ V {2u a - p u a - a - -(u a - a ) 2 }, 
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u lp ^ 
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(B4) 



We also need the RG equations for the umklapp cou- 
plings themselves, which are, 



( a P y}P + iL^y 1 ?) 



(B5) 
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where pg. = /«*-/« and gg = 2%cg + (2/g - /« - /ft), 
a = p, cr. 



APPENDIX C: INITIAL VALUES OF THE 
COUPLINGS 

Upon changing to the band basis, the on-site Hubbard 
repulsion is transformed into a set of interactions between 
the different bands. Using Eq. |2.3|, we have 



H int = -U^2 : c| T (x)c lT (x)c^(a;)c 4i (a;) 

i 



ijkl 



where A ijkl = ^S^S^S^Sr, 



(CI) 



(C2) 



Consider first OBCs. After linearizing around the 
Fermi points, each operator is split into left and right 
moving pieces. In terms of these, 



Hint =-U^2Y^ A ^ kl : ^ki^P2j^kkl^ 
ijkl P l= ± 
p i(-Plk F1 +P2k F2 -P 3 k F3 +P i kF i ) 



PM 



(C3) 



Now we are ready to compare the coefficients of inter- 
actions in Eqs. 2.29, C3. For example, comparing the 
coefficient in front of the term ipR^ip ' mi^uity 'li] yields 
the relation 



2 JlJ 



UA 



UA 



1133 ■ 



(C4) 



Comparing the coefficients of V'ffijV'.RijV'lij.V'Lij. gives 



(C5) 



We can then solve for the initial values of the forward 
couplings. 



where B. tj = ^ \S ml \ 2 \S m j\ 



A straightforward computation gives 
\ „ 1 . 



N + l 



(1 



Ji+3,N+1 



(C6) 



(C7) 



By similar comparisons, we obtain the initial values of 
the Cooper couplings, 



WAi 



2UBi 



(C8) 



where in the last step we use the fact that the transfor- 
mation matrix Sij in Eq. 2.5 is real for OBCs. 

For PBCs, similar results can be obtained by this 
method. Taking care to note the different conventions 
for left and right movers for PBCs (see Eq. 2.17), the 



(B8) initial values of the forward and Cooper couplings are 



f? 



4/£ 



2UAnjj, 



c^j 4c^ 2U Aijij , 



(C9) 
(CIO) 



where i = —i. From Eq. 2.4, one can compute these 
initial values. These are 

2U 
~N 



2U 



N 



L y ^13 



(Cll) 
(C12) 



If the number of chains is even, we also need the initial 
values of the transverse umklapp couplings. In fact, the 
initial values are the same as those for the forward and 
Cooper couplings: 

2U 



4u„ lp 



gi? 

y 13 jy ' 



ear = K - 



2U 
~~N~' 



(C13) 
(C14) 



Note that, in all case, the initial values of the rescaled 
couplings in Eqs. [T7-3.10 (without the tildes) are ob- 
tained by multiplying the factor l/7r(uj + Vj)- 
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APPENDIX D: REPRESENTATIONS FOR THE 
KLEIN FACTORS 



The Klein factors defined in Eq. 4.9 satisfy the com- 
mutation relations, 



{Via, Wjp} = 2 <5ij(5a/3. 



(Dl) 



In order to bosonize the relevant interactions in Eq. |4. 18[ , 
we need to prove that the products of the Klein factors 
in different terms commute with each other. Then, they 
can be simultaneously diagonalized with a specific choice 
of representation. 

T he pr oducts of the Klein factors for the first term in 
Eq. f4.18| are, 



VdiVdiVdiVdi = VdiVdTVdiVdi = 1. (D2) 
For the second term, the Klein factors we need are, 



ValValVblVbl = ValVa^Vbl = 9- 



(D3) 



A simple computation gives g 2 = 1. Thus, all the prod- 
ucts of Klein factors in the above equations commute 
with each other. It is therefore consistent to choose the 
trivial representation g = 1. For the cases of interest, 
then, no special signs or auxiliary Fermion fields are nec- 
essary in the bosonized Hamiltonian. 
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